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A METHOD FOR CALCULATING 
EXTERNALLY BLOWN FLAP NOISE 


Martin R. Fink 

United Technologies Research Center 


SUMMARY 


A method is described for calculating externally blown flap noise as a 
sum of several simple basic noise components. These components are (l) com- 
pact lift dipoles associated with the wing and flaps, (2) trailing edge noise 
associated with the last trailing edge, and (3) quadrupole noise associated 
with the undeflected exhaust Jet, deflection by a flap surface or nozzle de- 
flecting surface, and the free jet located downstream of the trailing edge. 
These noise components are combined to allow prediction of directivity and 
spectra for under -the -wing (UTW) slotted flaps with conventional or mixer 
nozzles, UTW slotless flaps, upper -surface -blowing (USB) slotless flaps, and 
engine-in-front -of -the-wing slotted flaps. The development of this method 
as part of a four-year effort under this Contract is described. A digital 
computer program listing is given for this calculation method. 

Directivities and spectra calculated by this method, the current NASA 
ANOP method, and a method developed by Lockheed -Georgia Co. are compared 
with free -field data for UTW and USB configurations recently tested by NASA 
Lewis Research Center. These data had not been part of the data base used in 
development of these three methods. The UTRC method best predicted the 
details of the measured noise emission, but the AJTOP method best estimated 
the noise levels directly below these configurations. 
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INTRODUCTION 


Noise generated by stationary solid bodies in the presence of the 
turbulent airflow in fan ducts is a major part of the noise of installed 
turbofan aircraft engines. For example, acoustically treated splitters 
within the engine inlet and exhaust ducts can attenuate turbomachinery noise 
but produce noise at their edges. Internal struts, necessary for structural 
support of the engine and splitters, are likely to be immersed in high-veloc- 
ity turbulent engine airflows. Turbofan stator blades are subjected to 
fluctuating wakes produced by the fan rotor blade. For STOL aircraft, 
externally blown flaps deflect engine exit airflow in order to generate 
additional wing lift force at low flight speeds. In all these cases, a solid 
surface of finite extent is scrubbed by airflow containing velocity and 
pressure fluctuations. The same basic aeroacoustic mechanisms should be pre- 
sent for all of these examples . 

The subject of this report is prediction of noise caused by externally 
blown flaps. This comprised the major portion of the investigation conducted 
under this four -year contract. Experiments with simple configurations were 
utilized to develop simple analytical models of several noninteracting noise 
components. Some components could, in concept, be calculated by rigorous 
methods. Such a theory would be extremely complex and is not presently 
available for externally blown flaps. Therefore a semi-empirical noise 
component method for externally blown flap (EBF) noise was developed instead. 
This UTRC method can be applied to under -the -wing slotless or slotted flaps 
and conventional or mixer nozzles, upper -surface -blowing configurations with 
uniform exhaust' flow, and engine -in-front -of -the -wing installations. Other 
prediction methods also exist for some of these types of EBF configurations. 
All of these methods are described and are evaluated herein by comparing their 
predictions with NASA -designated free field data for a range of EBF configu- 
rations. These recent data had not been used in development of these methods, 
so that the comparison -would not be biased. Crosscorrelation measurements 
to identify noise source locations are reported in Appendix D herein. 

A companion final report describes the remainder of the contract effort, 
devoted toward the development of a method for predicting and reducing noise 
radiation from gas turbine engine struts and splitters. 

Editorial review of this Contractor Report was performed at NASA Lewis 
Research Center and by colleagues within UTRC in order to assiire clarity of 
ideas expressed. Drafts of this report were also reviewed by persons at 
NASA Lewis Research Center and at Lockheed-Georgia Co. who have developed 
competing methods , in order to assure correct evaluation of their EBF noise 
prediction methods. 
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DESCRIPTION OE TIlE UTRC EBF NOISE PREDICTION METHOD 


Development of Analysis Under This Contract 

The work conducted under this Contract was based on an earlier noise 
component method developed at UTRC and first presented in reference 1. That 
method had regarded under -the -wing (UTW) externaly flown flap (EBF) noise as 
a sum of three types of noise components. Scrubbing noise as described 
therein was taken to be a lift dipole noise acting on the wing and flap 
panels scrubbed by the exhaust Jet. The noise generating process was taken 
as that described by Sharland (reference 2) for noise produced by a turbulent 
boundary layer on an airfoil surface. Such noise was assumed proportional to 
velocity raised to the sixth power, surface area scrubbed by the Jet, and 
surface pressure fluctuation squared. Pressure fluctuations on the flap 
surfaces were assumed equal to those for the free-Jet mixing region, an order 
of magnitude larger than those for an attached boundary layer. It was shown 
that equation (lO) of reference 2 predicted lobes of flap -generated lift 
dipole noise that generally matched the available EBF noise data at low ex- 
haust velocities for directions above and below the flaps. Another noise 
component was quadrupole noise generated by deflection of the exhaust Jet. 
This noise contribution was evaluated from the data of reference 3 for a Jet 
deflected by a large flat surface. Slotted and slotless UW configurations 
were noted to have approximately equal peaks of quadrupole -like noise at 
shallow angles above and below the deflected Jet. Measured deflected Jet 
noise radiation was therefore assumed to apply at both these directions de- 
spite lack of an explanation for the upward -radiated noise. The sum of these 
quadrupole and lift dipole components generally matched the available data 
except for underprediction of noise measured in the upper forward quadrant. 
The only noise generating process that is strongest at directions opposite to 
the deflected flap trailing edge would be trailing edge noise. This addi- 
tional noise component therefore had been assumed in order to match measured 
UTW directivity shapes in all quadrants . 

The first fiscal year's work under this Contract, reported in references 
4, 5j and 6, consisted of experiments directed toward evaluating the several 
available theories for incidence fluctuation noise, trailing edge noise, and 
scrubbing noise. Results for incidence fluctuation noise were presented in 
reference 5. That portion of the Contract effort, and additional Contract 
effort directed to predicting noise generated by struts and splitter rings 
within turbofan exit ducts, is summarized in reference 7. The study of 
trailing edge noise was also presented in reference 5* Those data, and 2/3 
of the summary data given by Hayden in reference 8, were found to substan- 
tiate the functional dependence of trailing edge noise developed by 
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Ffowcs Williams and Hall (reference 9) and by Chase (reference lO) . They did 
not agree with the dependence stated by Hayden in references 8 and 11. 

(Hayden has published a rebuttal in reference 12 which presents his data for 
the characteristic decay flow regime, which agrees with his equation and was 
not shown in reference 4. That comparison did not contain his data of 
reference 8 for a large range of velocities in the potential core or radial 
decay regimes, shown in figures 19 and 20 of reference 4 and in figure 6 of 
reference 5 • ) 

The study of scrubbing noise, published in references 4 and 6, had a 
major effect on the EBF noise calculation method developed under this 
Contract. It was shown by crosscorrelation experiments that regions of the 
flap surface which experienced large fluctuations of static pressure were not 
strong generators of noise. This result is frequently found in crosscorrela- 
tion experiments, as noted in Siddon's basic studies (reference 13) of sur- 
face-radiated noise. It contradicted the major noise process assumed in 
reference 1. Variation of maximum normalized crosscorrelation coefficient 
with chordwise position and variation, with chordwise position, of the times 
at which maximum fLrst derivative of crosscorrelation coefficient occurred, 
were utilized to develop a different description of this noise process during 
the second fiscal year. This concept also utilized the then -recent concepts 
of large-scale vortex structure in a UTW exhaust Jet as influenced by feed- 
back from the downstream solid surface (reference l4) . These vortexes, con- 
vected past the wing and flap surfaces, were assumed to induce lift force 
fluctuations on those surfaces. Strength of the lift force fluctuation and 
therefore the noise radiation depended on distance between the vortex trajec- 
tory and the surface. The first year's study of incidence fluctuation noise 
provided analytical Justification of the empirically assumed asymptotic 
spectrum slopes at large and small Strouhal numbers for this noise concept. 
The explanation developed at MSA (reference 15) to explain the measured 
variation of UTW noise with exhaust velocity to the 6.7 power rather than 
sixth power was also included for this component. That is, local exhaust 
velocity as calculated for the nozzle exit Mach number and flap trailing edge 
location was taken as the relevant flow velocity. 

During this second fiscal year, a method was developed for predicting 
the increased quadrupole noise caused by an exhaust Jet impinging against a 
flap surface. This calculation procedure used the MSA -developed method of 
reference l6 for noise radiated by an isolated Jet. It was shown that if 
this noise was increased by an amount proportional to sine squared of deflec- 
tion angle, and the directivity pattern was rotated through that angle, the 
data of reference 3 could be closely matched. OASPL directivities and I /3 
octave spectrum shapes were then calculated for most of the small-scale EBF 
configurations tested at MSA Lewis Research Center. A detailed description 
of the calculation method, and comparisons between predictions and data, were 
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given in reference 17 . Sketches of the twelve small-scale EBF ' 
which form the data base from which this UTRC method was develoi 
in Table I. 


TABLE I 

EXTERNALLY BLCWN FLAP MODELS FOR COMPARISON 
OF MEASURED AND PREDICTED NOISE 






UTW, 60° flap 
slotless wing 





UTW, 60° flap 
mixer nozzle 








During the third fiscal year, the calculation method was modifed to 
account for the fact that all the reference data base had been obtained under 
nonfree-field test conditions. Crosscorrelations among surface pressures on 
flap upper and lower surfaces, and far field acoustic pressures, were conduc- 
ted for both UTW and upper surface blowing (USB) configurations. These 
results generally validated the assumed conceptual model. Noise calculations 
were conducted for large-scale UIW and USB EBF configurations tested at NASA 
Lewis Research Center and corrected to free field. Comparisons with these 
UTW and USB data, for which the exhaust jet was either a nominally half- 
scale cold jet or a TF-34 engine, were given in reference l8. The UTRC 
method for calculating EBF noise was frozen at this time, except for two 
changes. Forward flight effects were included later this third year. Pre- 
diction of noise radiation from UIW slotless wings was modified during the 
fourth year because of poor agreement found during the third year. At this 
time the calculation procedure consisted of two separate time-sharing digital 
computer programs, one for OASPL directivity of each noise component and of 
their sum, and one for spectrum. Considerable manual effort was needed in 
preparing inputs for each program. Forward flight effects on local flow 
properties, surface pressure fluctuations, and far field acoustic pressures 
were measured in an acoustic wind tunnel for a range of exhaust velocities 
at different ratios of flight velocity to exhaust velocity. A method was 
developed for calculating the resulting effects of forward flight on various 
noise components. One unexpected result was that the spectrum of USB sur- 
face-radiated noise is decreased in amplitude and shifted to higher fre- 
quency. Full scale spectra, at high-annoyance frequencies, also decrease in 
amplitude as frequency is increased. The result is a large predicted noise 
reduction at low frequencies, where the contribution to annoyance is small, 
and negligible forward flight effect on high -annoyance noise. This predicted 
behavior is substantiated by data. The UTRC method for predicting EBF noise 
was modified to include these predicted effects of forward flight. Results 
obtained during the third fiscal year were presented in reference l8. Papers 
containing portions of the results were given in references 19 and 20. 

Activities conducted on EBF noise during the fourth (final) fiscal year 
included a test program to examine USB noise source locations by crosscorre- 
lation of local velocity fluctuations and far-field acoustic pressures. The 
noise calculation method developed by Tam and Reddy in reference 21 had 
assumed that one of the dominant components of USB noise is the flow mixing 
process in the highly sheared region downstream of the trailing edge and below 
the deflected exhaust jet. Instead, the crosscorrelation results presented 
in APPENDIX D herein Indicated that the noise was associated with large- 
scale structure of the exhaust jet that existed in the wing upper-surface 
boundary layer. These convected eddies produced noise only as they moved 
past the trailing edge. As another task, the several parts of the computer 
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programs for EBF noise as used in late 1975 were combined into one FORTRAN 
program, and a geometric error in calculating the effect of azimuth angle was 
corrected. This program calculates OASPL and I /3 octave spectrum for ranges 
of polar angle and azimuth angle specified by the input, for UTW configura- 
tions that can have conventional or mixer nozzles and slotted or slotless 
wings, USB configurations, or engine -in-front -of -wing configurations. Calcu- 
lations were compared with free -field data obtained by NASA Lewis Research 
Center and with limited additional data. Results for this fourth year are 
presented herein; portions of the results were presented in references 19 and 
22 . 


Assumptions Within UTRC Method 


Basic Concepts 

The noise component method described in this report calculates the total 
noise as a sum of several components which are acoustically but not aerody- 
namically independent. In its initial development, each component had the 
analytically predicted functional properties of noise radiation from a simpli- 
fied physical situation. If a fundamental analysis existed for noise radia- 
tion from surfaces in nonuniform flow, and if spatial distributions of mean 
velocity and turbulence spectrum were known for each EBF configuration, then 
the noise radiated by each component could be calculated rigorously. This 
information is not presently available. Empirical constants therefore have 
been developed in the analytical description of each noise component. These 
constants represent a combination of both the unknown local exhaust proper- 
ties of the deflected Jet (mean velocities, turbulence rras amplitudes, and 
turbulence integral scale lengths) and the unknown interactions between 
adjacent flap segments. From this viewpoint, the same basic simple noise 
components are assumed to occur for UTW slotted or slotless EBF configura- 
tions with conventional or mixer nozzles and for USB configurations. Ampli- 
tude of each noise component is calculated for the specific geometry of each 
case, and the resulting OASPL 's and spectra for all components are summed as 
statistically independent quantities. 

Noise radiation from surfaces in turbulent flow can have different 
behavior depending on the relative sizes of the chord, turbulence, and acous- 
tic wavelength. One convenient approximation used in reference 25 was to 
compare measured noise radiation patterns with those predicted for two 
limiting cases: very small and very large chord relabive to turbulence scale 

length and acoustic wavelength. The EBF noise prediction method given hereir^ 
denoted as the UTRC method, approximates the actual noise directivity pattern 
of surface -radiated noise as a sum of those two limiting cases. The limit of 
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very small chord is called fluctuating lift noise, and that of very large 
chord is called trailing edge noise. Functional dependence of noise ampli- 
tude is taken from theory; absolute level had been picked to match selected 
EBF data. The empirical agreement v?ith data was improved if another compo- 
nent, having the directivity shape and general behavior of jet mixing noise 
rotated through the exhaust deflection angle, was also included. 

Typical directivity patterns and relative amplitudes are sketched in 
figure 1 for each of these assumed noise mechanisms as it occurs for UW and 
USB configurations. Fluctuating lift noise, also called scrubbing noise or 
inflow noise, is sketched in the upper part of this figure. It is defined as 
an acoustically compact lift dipole noise oriented perpendicular to each 
chordwise segment of the wing. As shown experimentally in references 4 and 
6, local fluctuations of airfoil loading occurred which were coherent along 
the width of the region scrubbed by the exhaust jet. They moved downstream 
along the surface at the eddy convection velocity of about 80 ^ of local maxi- 
mum velocity. Local pressure fluctuations induced by this process resemble 
those caused by discrete vortices convected past an isolated airfoil. Ampli- 
tudes of these pressure fluctuations are small compared with pressure fluc- 
tuations generated by the shear -layer mixing process and impressed onto the 
ad.iacent surface. Thus the regions having strongest surface pressure fluc- 
tuations generally do not have strongest local noise source strength. This 
result that local dipole noise radiation is not necessarily proportional to 
local surface pressure fluctuation has been found in other studies (e.g., 
reference 13) of noise generation processes investigated by crosscorrelation 
techniques. 

As sketched in the upper part of this figure, far -field lift fluctuation 
noise from the undeflected portion of a wing and from each separately deflec- 
ted chordwise flap segment was represented by a separate dipole. Because an 
UTW flap (upper left sketch) is deflected into the jet exhaust and therefore 
closer to the hypothetical outer edge of the jet, the dipole associated 
with the aft flap segment is relatively strong. In contrast, an USB flap 
(upper right sketch) is deflected away from the jet exhaust. Its strongest 
assumed dipole tends to be that from the undeflected part of the wing. 

Trailing edge noise, sketched in the second row from the top, has a 
directivity pattern that is strongest directly upstream from the deflected 
trailing edge. Other properties of trailing edge noise are discussed in ref- 
erences 4 and 5- Noise generated at intermediate trailing edges such as that 
of the undeflected forward part of the wing was neglected. Thus the calcula- 
tion method developed here does not predict more trailing edge noise for mul- 
tiple slotted flaps than for single slotted or unslotted flaps. The trailing 
edge noise component was included because its directivity pattern tends to 
fill the gap in the upper forward quadrant between lobes of fluctuating lift 
noise, and thus produce better agreement between predictions and data. 
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Quadrupole noise, sketched in the third row, is represented as a sum of 
three components for both UTW and USB installations. One component is the 
jet mixing noise from undistorted, undeflected parts of the exhaust jet. 

This noise is calculated for the isolated exhaust nozzle and increased 2 dB 
to account for reflection of noise by the wing surface. Such noise is 
radiated beneath the UTW and above an USB configuration. UTW installations 
also generate a quadrupole noise from the region where the jet is deflected 
by the flap. This noise was called impact noise in references 24 and 25 . For 
conventional UTW installations, this quadrupole noise generated by deflection 
of the jet is stronger than that from the undeflected jet. This noise is 
radiated both above and below slotted flaps. USB configurations generate a 
similar increase of noise, radiated only above the wing and .flap, caused by 
deflection of the exhaust jet through the nozzle roof angle or cant angle. 

Both types of configurations produce additional quadrupole noise from the 
shear layer that forms beneath the deflected exhaust jet downstream of the 
trailing edge. Finally, as sketched in the lower right portion of the figure, 
external flow deflection devices which produce attached flow of the USB ex- 
haust jet to the deflected flap can radiate noise. Flow deflectors can gen- 
erate high-frequency dipole noise which is usually shielded by the wing, but 
some of this noise can be radiated below the wing. 

This empirical analysis models the exhaust jet as a line of discrete 
vortexes at the jet edge. Discrete vortices convected along an airfoil are 
known to induce local loadings concentrated near the vortex. The loading 
strength is a function of vortex chordwise position and varies approximately 
inversely with distance between the vortex and surface at constant chord. 

Too close a spacirg will cause viscous dissipation of the vortex, reducing 
the scrubbing noise. If the spectrum of vortex strength is that for jet tur- 
bulence and the lift force response is that for a discrete vortex in subsonic 
con^jressible flow, power spectral density of an acoustically compact source 
should vary as frequency squared at low reduced frequencies and frequency to 
the - 7/3 power at high reduced frequencies. The resulting l/3 octave slopes 
of 9 dB/octave and -4 dB/octave for low and high reduced frequencies, along 
with dipole directivity and dependence on local velocity to the sixth power, 
are typical properties of observed fluctuating lift noise. 

The concept of fluctuating lift noise as a lift dipole noise radiated on 
both sides of a wing and deflected flap, but generated by hypothetical coher- 
ent vortices in the exhaust jet on one side of the wing and flap, is funda- 
mental to the prediction method. This explanation was validated by tests 
described in reference 4 in which far-field spectra measured on both sides of 
a wing were compared. If this concept is correct, spectra measured on the 
side opposite from the jet should have a sum of trailing-edge noise and lift- 
dipole noise. Spectra measured at the same angle from the wing chord plane. 
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but on the side adjacent to the exhaust jet, should be the sum of those two 
noise processes plus jet mixing noise radiated directly to the far field and 
jet mixing noise reflected from the wing to the far field. To test this 
assumption, spectra measured both above and below an undeflected wing tangent 
to an exhaust jet were compared in figure 3^ of reference 4. Spectra for the 
isolated nozzle at the same direction angles and pressure ratios also were 
shown. Part of this figure for 250 m/sec exhaust velocity is reproduced as 

figure 2 herein. At directions 60° and 90° from upstream, spectra measured 

on the side adjacent to the jet (dotted line) were closely predicted 

(circles) by taking the spectra for the nozzle alone (dash line), adding 3 dB 

for reflection of jet noise from the wing, and adding that sum to the spec- 
trum measured on the shielded side (solid line). Further comparisons of 
spectra measured on the shielded side at different directions and exhaust 
velocities, described in reference 4, substantiated that those spectra were 
a sum of two simpler spectra. One of these components had maximum amplitude 
at a relatively low frequency, decayed rapidly in amplitude at higher fre- 
quencies, and had the fifth-power velocity dependence and cartioid direction- 
angle dependence of trailing edge noise. The other component had a broader 
spectrum shape with less rapid high-frequency decay. Its amplitude varied 
with velocity to the sixth power, and its directivity was that of a lift 
dipole. This latter noise component is what has been described herein as 
fluctuating lift noise. 

UTW Slotted Wing 

Geometric properties of a UTW slotted wing are given (figure 43) by the 
number of flap slots, leading edge ordinates for the wing and each flap panel 
relative to the center of the nozzle exit, deflection of the wing chord line 
relative to the nozzle centerline and of each flap relative to the wing chord 
line, and chord of the last flap segment. If aerodynamic data are available 
for the effective turning angle of the flap segment, deflection angle of the 
last flap segment is taken as this aerodynamic angle rather than the geom- 
metric angle. Chord of the wing, and of all but the last flap segment, is 
taken as the distance from its leading edge to the next leading edge. How- 
ever, trailing edge location is calculated from chord and the input deflec- 
tion angle. The computer program can accept up to four chordwise flap seg- 
ments (quadrupole slotted flaps). UIW wings with retracted flaps are dis- 
cussed in the following section entitled "UTW Slotless Wings". 

A hypothetical line of vortices is assumed to induce lift force fluc- 
tuations on the wing and flap segments. Vortex trajectory is taken as a 
straight line, parallel to the nozzle centerline, extending downstream from 
the nozzle lower lip until it gets within half a diameter of the flap surface. 
If the wing or flap extends below the nozzle centerline, the vortex 
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trajectory becomes parallel to the flap chord line and displaced half a 
diameter away from it . 

Lift force flactuation caused by a vortex convected past an airfoil is 
known to vary inversely with distance of closest approach. Noise radiation 
is assumed to vary with force fluctuation squared and therefore directly with 
the product of nozzle diameter and total chord and inversely with average 
distance squared. For each chordwise segment, the average value of inverse 
distance squared, h”^, is taken as the average of this quantity evaluated at 
the segment leading and trailing edge. Spreading of the exhaust jet by the 
deflected flaps was neglected, so the mean square acoustic pressure from each 
segment is assumed proportional to the product of segment chord and nozzle 
diameter divided by far field distance squared. Noise amplitude caused by a 
vortex distance of one diameter was selected to match data from reference 4 
for a wing with retracted flaps. By varying the assumed trajectory, the vor- 
tex minimum displacement of half a nozzle diameter was inferred as that value 
which would predict the measured noise increase reported in reference 26 in 
changing the flaps -retracted to the approach configuration for a doable 
slotted flap. (Absolute amplitudes for noise radiation were later decreased 
because the data given in reference 26 are not free -field.) Thus for the 
wing and for each of n flap segments 

Kn= 1 Xto-'^{Cn/c) (hn/c)'^ (1) 

and the total overall fluctuating lift noise is given by 

2 

OASPLl= loiog (cD/r^) (U,/vf ( 2 ) 

As with the basic concepts of the NASA ANOP method of reference 27 for 
EBF, fluctuating lift noise was assumed to vary with impingement velocity 
raised to the sixth power. Impingement velocity was defined as the maximum 
velocity in an isolated axisymmetric exhaust jet at the axial distance of the 
impingement point. For a flap configuration that extends below the nozzle 
centerline, the impingement point is the position where the centerline inter- 
sects the flap. For flaps that are not deflected that far into the exhaust 
jet, the impingement point is taken as the trailing edge of the Last flap 
segment. The ratios of impingement velocity and trailing edge velocity to 
nozzle exhaust velocity were calculated from the NASA -developed equation ( 3 ) 
of reference 15- 

Uj/ V = {1 + [o.l4 (X/D) (1 + ( 3 ) 


li 


Large ratios of flap length to nozzle diameter have been found to generate 
less noise than this procedure would predict. Lift fluctuation noise there- 
fore was arbitrarily multiplied by the ratio of trailing edge velocity to 
irapingeirient velocity squared. Equivalent nozzle velocity for sin unmixed 
coajiial jet was calculated from equation (5) of reference 27» the NASA ANOF 
method for EBF. 

Trailing edge noise is proportional to the product of turbulence 
integral scale length and spanwise distance along the edge. Both of these 
dimensions are proportional to nozzle diameter. Amplitude of trailing edge 
noise is proportional to the ratio of diameter squared to far-field radius 
squared, and to trailing edge velocity raised to the fifth power. Directiv- 
ity is given by cosine squared of half the angle from the last flap segment's 
upstream direction. 

OASPL^ = 10 log 0.5X10'® {/U|®/ 0 P,ef) (D/rf COS^'^ C0s2[ (0 + S f /Z)] (1. ) 

Normalized 1/3 octave spectra taken from equations (ll) and (12) of 
reference 17 are used for fluctuating lift noise and trailing edge noise. 
These spectra have analytically justified asymptotic slopes at small and 
large Strouhal numbers, but nevertheless they are empirical curves based 
on published normalized spectra. The sole justification for use of these 
equations is the good agreement with measured normalized EBF spectra given 
in reference 1. 


(SPL|/ 3 - 0ASPL)l= 10 log 0.037 St‘^(Sl®^®+0,008)“^ (5) 

{SPL|/ 3 - OASPL)^- 10 log 0.029 (6) 


Calculation of the several kinds of quadrupole noise is based on the 
NASA AJ'JOP method of reference l6 for quadrupole noise of an isolated .axis^nn- 
raetric subsonic exhaust jet. OASPL is calculated by an explicit equation, 
with an empirical correction for refraction at direction angles close to the 
downstream centerline. Normal to the jet. 


OASPL 


90 


= 141 +10 log [m/'^ (i + 0.010 m/®)"' 
10 [3m/-®( 0.60 + Mj®-®)‘'-l] log (Pj 


\ + 10 
/Pq) 


log [(A/r^) (Po/p,s/(ao/a,sA)^] + 

(7) 


For direction angles, relative to the deflected Jet upstreaiii direction, 

smaller than M xl 50 °, 

(.] 


OASPLj OASPL gQ “30 log j^l + Mj; (l + Cos 


(B) 


where the convective Mach nu.-!iber M was taken as 0.62M . 

c J 


For relative direction angles between Mjxl50^ a-iid l80 , the negative quantity 

AOASPj = - 18 [l-(6-MjS/30)] 

is added to OASPLj calculated from equation (8) to account for refraction 
within the jet. Deflection of the exhaust jet by a flap segment is assumed 
to add a quadrupole noise term given by 6 sine squared of the effective turn- 
ing angle (input as the last flap angle) multiplied by noise of the isolated 
jet. The sum of these 'tv; o quadrupole terras 

OA5PLo= OASPLj+ 10 log [( 1+6 sin^S f)(> + COS^ (^)/2] (lO) 

has its directivity defined relative to the deflected -flow centerline. 
Quadrupole noise is also assumed to be generated by the portion of the ex- 
haust jet downstream of the last flap segment trailing edge. This noise is 
arbitrarily taken as that for an isolated jet having a diameter equal to the 
nozzle diameter and exhaust velocity equal to the calculated maximum velocity 
at the trailing edge. It is then multiplied by the same factor, 1+6 sin^ 
(deflection angle), used with the jet deflection process. This quadrupole 
noise component is also referenced to the deflected-flow centerline. All of 
these quadrupole noise components radiate equal intensities above and below 
the deflected jet. This assumption was made because experimentally deter- 
mined UTW directivity patterns contain equal-amplitude peaks of quadrupole 
noise above and below the deflected jet. 


Quadrupole noise at directions below the jet, termed impact noise in 
reference 2^4, was taken in that study as equal to the noise measured at the 
same angle beneath a jet deflected by a large solid surface. Acoustic re- 
flection from that simplified configuration raises the noise caused by the 
undistorted portion of the jet between the nozzle exit and the impingement 
region, and by jet distortion and deflection in the impingement region. 
Whether slotted flaps were assumed to transmit or reflect this noise, it was 
necessary to postulate an additional quadrupole noise process having approxi- 
mately the same strength and orientation. Their sum would match UTW data at 
directions 20° to 4o*^ above and below the deflected jet. The quadrupole 
noise arbitrarily attributed to the exhaust jet downstream of the trailing 
edge was hypothesized for this reason. 


UTW quadrupole noise is important at direction angles greater than 90° 
from the upstream direction of the d.eflected jet. V/ithin this region, nor- 
malized spectrum shape for an isolated jet changes greatly (figure 5 of ref- 
erence l6) with direction angle. Spectrum shapes for directions near peak 
OASPL amplitude have a more rapid high-frequency decay tlmn that for more 
upstream directions. It was found that spectrum shapes which were correct 
for an isolated jet at peak OASPL underpredicted the measured high-frequency 
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noise for UIW slotted flaps at approach flap position. Closer agreement was 
obtained by use of the normalized spectrum for an isolated jet at l40° direc- 
tion angle for all angles between l40° and l8o°, 

(spl,/ 3 - OASPL)q = 10 log 0.1 (st''^^'^+0.ll)''^ (11) 

The contribution of quadrupole noise to total UTV/ noise below the wing 
decreases as direction angle is decreased from this range. Rather than in- 
clude the complexity of a spectrum shape that varies with direction angle, 
over a direction range where this contribution becomes small, this one nor- 
malized spectrum for quadrupole noise is used at all directions. 

As part of the development of the UTRC method, calculated variations of 
UTW noise with changes in configuration geometry were compared with available 
data. Measured effects of axial and vertical displacement of a double 
slotted flap approach configuration (reference 26 ) were shown in reference 17 
to be closely predicted. Measured effects of doubling and halving the nozzle 
diameter of that model, at constant axial and vertical position of the nozzle 
upper edge, were also (reference l8) correctly given. Although based on data 
for that double slotted flap at 20 ^ and 60 ° flap deflections, this fJTW noise 
prediction method was shown in reference l8 to predict closely the measured 
noise from a 4o° deflection triple slotted flap. 

Forward flight effects on UTW fluctuating lift noise and trailing edge 
noise, taken from references l8 and 20, are represented as a decreased ampli- 
tude given by the ratio of relative velocity to exhaust velocity squared. 
Forvjard flight effects on quadrupole noise from the undeflected portion of 
the exhaust are assumed to vary with relative velocity ratio to the sixth 
power as with the NASA ANOF method of reference l6 for noise of an isolated 
jet. However, the increases of quadrupole noise attributed to impingement 
against the flap lower surface and to the jet dovmstream of the trailing edge 
are assumed to vary with relative velocity ratio squared. This reduced expo- 
nent was based on the data of reference l8. Relative velocity betv/een the 
exhaust jet and the surface against which it impinges is not affected by for- 
ward flight. These calculations provide the predicted effect of forward 
flight on EBF noise source strength in a coordinate system fixed to the air- 
frame. They do not include the effect of airframe motion relative to a 
ground-fixed observer. Corrections for this difference betv;een an aircraft 
flyover and an acoustic wind tunnel test are given in reference 27 . 
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UTW Mixer Nozzle 


Some UTW slotted flaps have heen operated with mixer nozzles to decrease 
the peak velocity and temperatiire at the deflected flaps. For this type of 
configuration, ratio of local velocity at the impingement point to exhaust 
velocity must he supplied as input. This quantity is arbitrarily defined as 
the average of the two largest local peak values of velocity ratio measured 
in the exhaust of the isolated mixer nozzle at the impingement -point axial 
location. The same velocity ratio is assumed to apply at the trailing edge. 
Nozzle diameter is taken as the hydraulic diameter of the mixer nozzle total 
exit area. Calculation of fluctuating lift noise and trailing edge noise 
then follows the method for UTW slotted flaps having conventional axisymmetric 
nozzles, except that the resulting levels are arbitrarily increased 3 dB. 

This increased amplitude, possibly caused by higher turbulence levels or 
larger-diameter core of high-velocity exhaust flow, was shown in references 
IT and l8 to be necessary for obtaining good agreement with data. 

This computer program does not calculate quadrupole noise for UTW mixer 
nozzles. Although OASPL amplitudes for a mixer nozzle generally match those 
for the unmixed exhaust jet, spectrum levels depend strongly on the mixer 
nozzle geometry. It is assumed that measured noise spectra and OASPL ampli- 
tudes are available for the isolated mixer nozzle at takeoff and approach 
flight speeds. This measured noise must be increased by 10 log (l+6 (veloc- 
ity ratio)® sine squared (flow deflection angle)) to account for deflection 
of the mixed exhaust jet. The resulting quadrupole noise should be rotated 
through the flow deflection angle, and all amplitudes at directions below the 
wing and flaps should be increased 3 dB. This increase had been stated in 
reference 17 and l8 to be caused by reflection of quadrupole noise from the 
wing and flap. However, improved agreement with data measured above the wing 
is obtained if upward -radiated noise is assumed to bo transmitted through the 
flap slots. The resulting quadrupole noise component (increased 3 dB below 
the wing and flap, not increased above them) should be added to surface- 
radiated noise given by the computer program. 

Upper Surface Blowing 

Geometry of a USB (figure 43) is idealized as two straight lines, one 
for the wing and one for the last flap. The wing is described by the upper 
surface coordinate at the nozzle exit plane and by the wing incidence rela- 
tive to the nozzle centerline. The flap is given by its trailing edge coor- 
dinates and by the aerodynamic turning relative to the wing chord line. 
Coordinates at the intersection of those two lines are computed as part of 
the solution. Effects of nozzle exit shape are neglected, and the nozzle is 
specified by its hydraulic diameter and its roof angle, cant angle, or exter- 
nal vane deflection angle relative to the nozzle centerline. Details of the 
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nozzle anH deflector shape have been found empirically to have significant 
effects on noise amplitudes; these can be estimated using the data given in 
reference 28. 

Trailing edge noise and fluctuating lift noise OASPL's are calculated in 
the same manner as for UIW slotted wings. However, amplitude of trailing 
edge noise was taken as twice as large for USB as for UTW. The vortex tra- 
jectory is assumed to be one hydraulic diameter above the wing. Lift fluc- 
tuation noise from the flap is arbitrarily taken as 1.5 times that for a 
hypothetical vortex trajectory one hydraulic diameter above the flap. Veloc- 
ity ratio at the trailing edge is calculated for a distance equal to the path 
length along the wing and flap. It was found that spectra calculated in this 
manner decayed less rapidly than the data for large Strouhal numbers and low 
exhaust velocities. Therefore the normalized spectrum shape for trailing 
edge noise was arbitrarily used with OASPL calculated from the sum of lift 
fluctuation noise and trailing edge noise. 

Quadrupole noise caused by impingement of the exhaust Jet against a 
deflecting surface is assumed to increase the quadrupole noise of an isolated 
exhaust Jet. This increase caused by roof angle, cant angle, or vane angle 
is given by the same factor as that for impingement of a UTW exhaust Jet 
against a deflected flap. The resulting quadrupole noise is assumed to be 
shielded by the wing and flap upper surface. It is rotated through the aero- 
dynamic turning angle and is calculated only for direction angles above the 
wing and deflected flap. Quadrupole noise from the portion of the exhaust Jet 
downstream of the trailing edge is increased by the same factor and rotated 
through the same angle, but it radiates to all directions. As with UTW, one 
normalized spectrum shape is used at all directions. It is likely that for 
direction angles at which USB quadrupole noise is important, improved agree- 
ment with data would be obtained if the normalized spectrum was assumed to 
vary with angle relative to the deflected Jet centerline as with an isolated 
Jet. Measured noise of each isolated nozzle and deflector, without the wing 
and flap, should be compared with quadrupole noise calculated for the deflec- 
ted Jet. Any increment of measured dipole noise associated with the presence 
of the deflector should be added to the predicted USB noise for directions 
above the wing. 

Forward flight effects on the sum of trailing edge noise and fluctuating 
lift noise was taken as a decrease of amplitude combined with an increase of 
frequency. As with UTW, the OASPL amplitude is taken proportional to rela- 
tive velocity ratio squared. However, normalized I/3 octave SPL is calcula- 
ted for Strouhal numbers that were multiplied by the quantity (l + flight 
velocity /exhaust velocity). This adjustment had been developed in references 
18 and 20. Forward flight effects on quadrupole noise are taken equal to 
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those for UTW. One reasonable validation of these predicted forward flight 
effects on USB noise -would he comparisons -with flyover data for the Boeing 
yC-l4 Advanced Medium STOL Transport, the only full-scale USB airplane, if 
noise from the installed engine proves to be sufficiently low. 

UTW Slotless Wing 

Geometry of a UTW slotless -wing is specified and calculated in the 
same manner as that for USB. Quadrupole noise and trailing edge noise are 
computed in the same manner as for a UTW slotted -wing. To obtain closer 
agreement between measured and predicted noise directivity, it was found 
necessary to add a noise component having the general shape and velocity 
dependence of a lift dipole oriented perpendicular to the deflected aft flap* 
However, increasing the flap chord was experimentally found not to increase 
this OASPL. This noise component was arbitrarily modeled as fluctuating lift 
noise with an amplitude three times that which would be calculated for the 
wing segment, rotated through the flap deflection angle. 

This special calculation of fluctuating lift noise is applied only if 
the slotless wing extends below the nozzle centerline. Any slotless UTW 
which does not meet this condition is regarded as a wing with retracted 
flaps. Geometry of a wing with retracted flaps is specified as that for a 
USB configuration having zero flap deflection. It could also be specified as 
a UTW single slotted flap having a flap with zero deflection and zero chord, 
located at the wing trailing edge. Both calculations will give the same 
result . 


Engine in Front of Wing 

This type of configuration has the wing chord line coincident with the 
exhaust nozzle centerline. The wing is conpletely immersed in the turbulent 
exhaust jet. For the one configuration- tested, maximum wing thickness was 
about half the nozzle exit diameter. Therefore the exhaust jet was substan- 
tially altered by the presence of the wing. Maximum local velocities at the 
trailing edge flaps probably were reduced by viscous decay along the wing. 
Spectra radiated by this model were characterized by a high-frequency decay 
rate larger than that for other EBF configurations . This rate of about 12 dB 
per octave for l/3 octave band spectra agreed with data for isolated airfoils 
in turbulent flow. 


A calculation method developed in references 4 and 5 for predicting 
noise radiated by isolated airfoils in three-dimensional compressible flow 
was applied to the undeflected wing panel and each flap segment. Each lift 
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dipole was assumed to act normal to the local chord. Calculated Jet mixing 
noise for an isolated undeflected exhaust Jet was added to the lift fluc- 
tuation noise. Width of the turbulent flow was taken equal to the nozzle 
diameter, and flow velocity was taken as the nozzle exit velocity. It was 
found that free-field spectrum could be matched if turbulence intensity was 
arbitrarily set equal to 7 percent and turbulence integral scale length was 
taken as one-eighth the nozzle radius. (The data had been measured with a 
reflecting ground surface, and a 10 percent turbulence level was cited in 
reference 17 to match those data.) Resulting calculated OASPL directivity 
tended to overpredict the data for the first and third quadrant at approach 
flap deflection. Closer agreement was obtained by neglecting the acoustic 
contribution of the last flap panel. It is possible that local flow velocity 
at that location was decreased by the presence of the wing within the Jet, 
causing a decrease of local noise radiation. 

Calculated noise radiation for this type of configuration should be 
applicable to interference noise of a wing-mounted propfan. Slipstream 
velocities for highly loaded propfans designed for cruise at Mach numbers 
near 0.8 would be comparable to exhaust velocities of lightly loaded high by- 
pass ratio turbofan engines. 


EVALUATION OF EBF NOISE PREDICTION METHODS 


Other EBF Noise Prediction Methods 

The method developed under this contract is evaluated in this section by 
comparing its predictions, and those from other openly published methods, with 
NASA-supplied data. These new data had not been part of the data base used 
in development of the EBF noise prediction methods compared herein. Capabili- 
ties and limitations of these methods are compared in Table II. 


18 



TABLE II - CAPABILITIES AMD LIMETATIONS OF 
EBF NOISE PREDICTION METHODS 
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tions having a uniform exhaust jet is currently limited to prediction of OASPL 
directivity; it does not predict spectra. One widely used method is the NASA 
ANOP (Aircraft Noise Prediction Program) method of reference 27 . This method 
is based on empirical curves drawn through data available in 1975; these data 
were corrected for ground reflection. Application to USB configurations 
having noncircular exhaust nozzles was not specifically described. However, 
noise from such configurations has been calculated by NASA using the ANOP 
method by replacing the actual nozzle with a circular nozzle of equal exit 
area. This same usage is made herein. 

Another method used for this evaluation is the GELAC (Lockheed -Georgia 
Co.) method developed for the FAA. It is a noise component method which con- 
tains separate calculations of surface -radiated noise and quadrupole noise 
from several geometric regions. The early version of this procedure (refer- 
ence 29) included the effects of many UIW geometric variables (flap vertical 
and axial position, nozzle pitch angle, and number of flap slots) not 
represented in the ANOP method. A more recent version (reference 30 ) has 
much less dependence on flap geometry and is more strongly directed toward 
USB configurations. Both versions exist as digital computer programs and as 
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hand calculation procedures. Predictions given herein were obtained by the 
hand calculation method of reference 29 for UTW and that of reference 30 for 
USB. Because the data base for the method of reference 29 consisted of mea- 
surements that were affected by ground reflection, predicted levels were de- 
creased 3 dB for closer agreement with MSA free field data. The earlier 
method of reference 29 represents all UTW noise as varying with exhaust 
velocity to the sixth power. This method was utilized because OASPL's and 
directivity patterns calculated from reference 30 were in very poor agreement 
with UTW data. It was later recognized that all quadrupole noise terms given 
in reference 30 were of the order of 10 dB below expected levels. Quadru- 
pole noise calculations developed in the GELAC method of reference 30 were 
scaled from levels predicted by the MSA ANOP method for jet exhaust noise 
(reference 16 ) . Jet exhaust noise for an isolated circular exhaust nozzle is 
given by equation ( 5 ) of reference 30, which agrees with the corresponding 
equation (6) of reference I 6 except for one item. The numerical value for 
the constant K used in reference 16 was l4l dB. However, this constant as 
specified on p. A -86 of reference 30 for hand calculations by the GELAC 
method was 134 dB. This same constant, 7 dB less than that specified for the 
MSA ANOP method for jet noise also appears within the GELAC computer program 
listing (line 110 of subroutine JET, p. A-24 of reference 30), Comparisons 
between calculated and measured OASPL were not shown within the GELAC docu- 
ment. Although not explained in reference 30, this change from the constant 
used in reference I 6 was included to provide closer agreement with unpublish- 
ed GELAC USB data and its effect on UTW predictions was not examined. 

A method was developed at NASA Lewis Research Center for calculating 
OASPL directivity of UTW configurations at directions under the wing in the 
flyover plane. Use of this method for slotless wings was first described in 
reference 24. The significant noise sources for slotless wings were assumed 
to be trailing edge noise and impact noise. Trailing edge noise was calcu- 
lated from local maximum velocity, boundary layer thickness at the position of 
maximum velocity, and width of the exhaust jet (all measured at the trailing 
edge). Impact noise caused by a Jet Impinging against a very large deflecting 
surface was obtained from data such as that of reference 3. This method was 
extended in reference 25 to include application to UTW slotted wings. An 
additional noise mechanism, inflow noise (called fluctuating lift noise in 
the UTRC method), is included for those configurations. This noise is assumed 
to result from the large-scale turbulence structure, present in Jet exhausts, 
moving past the flaps to cause fluctuations of lift force. The noise was cal- 
culated using local flow properties estimated at midchord of each flap panel. 
This NASA Lewis method currently does not predict noise spectra. OASPL 
directivities calculated for slotless wings by this method were provided by 
rir. D. J. McKinzie, Jr., of NASA Lewis Research Center. 
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II 


other calculation methods for USB had been developed by Hayden (refer- 
ence ll), Biller (reference 3l)j and Reddy (reference 21). Noise levels mea- 
sured directly beneath two USB models were compared in figures 10-12 of ref- 
erence 32 with predictions by the methods of references 11 and 31, and by an 
early version (reference 1?) of the method of this report. Upper surface 
length from the slot nozzle to the flap trailing edge differed by a factor of 
2 for these models. OASPL of both models was overpredicted roughly 10 dB by 
the method of reference 11 and underpredicted the same amount by the method 
of reference 31. For both of those methods, the predicted spectrum shape was 
more sharply peaked than the measured shape. Good agreement on level and 
spectrum shape was obtained with the early UTRC method. Those two prediction 
methods were not evaluated herein because they gave such poor agreement with 
those data. The more recent method of reference 21 was not available in 
sufficient detail to be applied. 


Designated EBF Model Configurations 

The EBF model configurations shown in figiires 3 and 4 were designated 
by NASA Lewis Research Center to be used for evaluation of noise prediction 
methods. These models were relatively small, with nozzle exit areas generally 
equal to that of a 5.2 cm (2 in.) diameter circular nozzle. All had been 
tested at NASA Lewis Research Center in an outdoor facility equipped with a 
mat of acoustic absorbing foam to give free-field data above 200 Hz. Data had 
been corrected for atmospheric attenuation. Comparisons were also made with 
nonfree-field data for specific unique configurations. 

Two designated models were under -the -wing configurations. Sketches of 
these models are given in figure 3 . These were a three-flap installation at 
both takeoff and approach flap deflection, for a range of sideline angles, 
and a slotless wing having the same lower surface contour as the three-flap 
approach flap shape, for only the flyover plane. Portions of these data had 
been presented in reference 19 . They were compared therein with predictions 
by the UTRC, NASA ANOP, and GELAC methods. 

The other three models, shown in figure 4, were upper-surface-blowing 
(over-the-wing ) configurations. One was the QCSEE configuration having an 
aspect ratio 2 slot nozzle, tested at nominal I/II .5 scale with l4 cm (5.5 
in.) equivalent nozzle diameter. Data were supplied for a range of flap length 
and sideline angle at takeoff flap deflection and a range of sideline angle at 
approach flap deflection. These data had been presented in reference 22 where 
they were compared with predictions by the UTRC method. Another was a I/I 8.5 
scale model of the TF-34 over-the-wing aspect ratio 4 slot-nozzle short-flap 
configuration of reference 33 with takeoff flap deflection. The third 
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configiiration (reference 3^) had a circular nozzle above the wing and an 
external jet-flow deflector to direct the exhaust jet toward the wing upper 
surface. This arrangement was expected to reduce adverse interference of the 
exhaust jet on the wing drag force d\aring cruise. 

Aerodynamic performance data (lift, thrust, and trailing edge velocity 
profile) at zero forward speed are available for all these configurations. 


Comparisons With Designated UTW Data 

Predicted and measured UTW noise radiation are compared in the following 
order: OASPL at the 90° polar angle 0° azimuth angle flyover position, 

general shape of OASPL directivity in the flyover plane, effect of azimuth 
angle on OASPL amplitudes, and normalized l/3 octave spectra. The discussion 
is directed primarily toward evaluation of the UTRC method, the subject of 
this report. 

Triple Slotted Flap, Approach 

This configuration had been tested at nominal 1/28 scale. Data were 
presented, and were compared with predictions by several methods, in refer- 
ence 19. The free-field spectra contained small peaks and valleys at non- 
uniformly spaced frequencies. These frequencies did not vary with exhaust 
velocity but were altered by changing the flap deflection. 

Calculated and measured OASPL directivities in the flyover plane for 
approach flap deflection and four exhaust velocities are compared in figure 5. 
At low exhaust velocities the UTRC and ANOP methods underestimated the mea- 
sured OASPL at 90° polar angle by 3 to 4 dB, and the GELAC method was about 

2 dB high. At the two higher exhaust velocities the UTRC and GELAC methods 
were within 1 to 2 dB and 1 to 3 dB, respectively; the AWOP method remained 

3 dB low. 


Both the ANOP and GELAC methods use a constant directivity shape that 
generally matched the data. In contrast, directivity shape calculated by 
the UTRC method varied with exhaust velocity. Except at polar angles near 
the deflected exhaust, the measured directivity shape did not change greatly 
with exhaust velocity. Measured OASPL amplitudes varied approximately with 
the 6.5 power of exhaust velocity. 
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Effects of sideline angle on measured and calculated OASPL for this 
approach flap configuration are shown in figure 6. These effects are plotted 
as changes in OASPL relative to that for the flyover plane, at two exhaust 
velocities, as a function of polar angle (angle relative to the nozzle cen- 
terline). The GELAC predictions (defined only for 90° polar angle) essenti- 
ally matcn the maximum noise reductions calculated hy the ANOP method. In 
contrast, the UTRC predicted maximum reductions are about twice as large and 
generally agree with the data. These calculations differ from those given in 
reference 19» which contained a trigonometric error. Details of the predic- 
ted shape do not precisely matcn tne data; maximum reduction is predicted near 
90° polar angle but occurred about 10° further downstream. For 60° sideline 
angle the predicted and measured maximum reductions of about 6 dB generally 
agree with the analytically expected dependence of surface -radiated noise on 
sine squared of tne sideline angle. For 85° sideline angle the measured and 
UTRC calculated sideline noise decrease roughly 10 dB to a floor set by 
quadrupole noise from the deflected Jet. Measured and UTRC calculated noise 
was increased several dB at 120° polar angle. This angle is within the re- 
fraction valley of quadrupole noise from the deflected exhaust Jet as viewed 
in the flyover plane, but is near peak amplitude of quadrupole noise when 
viewed from the side . 

Calculated and measured 1/3 octave spectra normalized with respect to 
OASPL, for approach flap deflection, are compared in figure 7« Frequencies 
are normalized as Strouhal number based on nozzle diameter and exhaust veloc- 
ity. Data are shown for 70° and 110° polar angle and 120 and 228 m/sec ex- 
haust velocity, in the flyover plane and at 85° sideline angle. For tne fly- 
over plane and large Strouhal numbers, the data points for 70° polar angle 
were about 5 d.B below those for 110° polar angle. OASPL was somewhat affec- 
ted by irregular peaks in the spectra measured at forward polar angles; 
actual 1/3 octave band SPL's for the two directions differed by about 10 dB 
near peak amplitude but were within 2 dB for Strouhal numbers above 2. All 
three prediction methods generally matched the normalized spectra for 110° 
and (not shown) 90° polar angles. Normalized spectra for the 85° sideline 
plane had about 8 dB data spread at large Strouhal numbers. The ANOP and 
GELAC methods bracketed tne data for small Strouhal numbers and generally 
matched the high-velocity spectra for large Strouhal numbers. Spectrum shapes 
calculated by the UTRC method varied because the relative amplitudes of dif- 
ferent noise components vary with exhaust velocity and polar angle . The 
resulting calculated spectra lie within a narrow envelope. The UTRC method 
was closer to an average of the data for large Strouhal numbers. It gave a 
more rapid spectrum roll-off at high frequency than did the other two methods. 
This stronger decay was in closer agreement with the data. 
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Triple Slotted Flap, Takeoff 


Calculated and measured OASPL directivities in the flyover plane for the 
takeoff configuration and four exhaust velocities are shown in figure 8. 
Calculation by the ANOP method used equation (l6) of reference 27, rather 
than equation (ll) of that reference, for predicting the variation of OASPL 
with flap deflection angle. Use of equation (l6) was recommended for config- 
urations in which the flaps extend relatively far into the high-velocity ex- 
haust Jet at small deflection angles. If the other equation had been used, 
calculated noise levels would have been 5*2 dB smaller. Measured amplitudes 
were underpredicted by all three methods. Near 90° polar angle the UTRC 
method and the ANOP method with the more favorable equation were about 5 dB 
low. The GELAC method of reference 29 was about 8 dB low. This ANOP equa- 
tion had been developed specifically from data for triple-slotted flap con- 
figurations with this type of geometry. In contrast with this poor agreement, 
these three methods were known to predict closely (within 2 dB) the OASPL 
directivity data of reference 26 for an UTW double-slotted flap at takeoff 
deflection. All three methods had been developed to match those data. That 
double -slotted flap configuration had a smaller total flap chord than the 
model triple-slotted flap, and about the same wing leading edge position and 
wing chord. At takeoff deflection, the' trailing edge of the double-slotted 
flap’s last flap segment was located above the nozzle centerline. From the 
viewpoint of the UTRC method, the deflected flap panels did not greatly dis- 
tort the exhaust jet so they did not produce much lift fluctuation noise. 
Nearly all the last panel of the triple -slotted flap extended below the 
nozzle centerline at takeoff deflection. None of these three methods correc- 
tly predicted the resulting increased noise levels. The GELAC method, which 
was the only method to overpredict measured levels near 90 ° polar angle for 
the approach configuration, gave the worst underprediction at takeoff. 

A fourth set of predicted curves, labeled NASA LEWIS, shows OASPL 
directivities calculated by D. J. McKinzie, Jr. of NASA Lewis Research 
Center by the method of reference 25. These calculations require extensive 
knowledge of local turbulence and mean velocity. Such calculations were 
supplied by NASA only for this takeoff configuration, where other methods 
gave worst agreement with data. These calculated amplitudes were within 2 dB 
of data near 90° polar angle. They generally matched the data from there to 
the deflected jet and were up to k dB high in the forward quadrant. This 
good agreement with data was achieved by use of estimated local flow properties 
evaluated from data for isolated jets, and illustrates the need for such flow- 
field information. As previously noted, this method has not yet been extended 
to prediction of spectra or of sideline directivity. 
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The effect of sideline angle on calculated and measured OASEL directiv- 
ity is shown in figure 9* Measured maximum reductions, relative to those in 
the flyover plane, were about 7 and 12 dB for 60° and 85 ° sideline angle, 
respectively. These reductions were more than twice those predicted by the 
ANOP method and were generally predicted by the UTRC method. The UTRC 
method also correctly predicted that the difference between OASPL in the 
sideline plane and the flyover plane would be positive at 150 ° polar angle . 
This polar angle is in the deflected exhaust jet's refraction valley for the 
flyover plane but not for the 60 ° or 85 ° sideline angles. 

Calculated and measured normalized l/3 octave spectra at takeoff flap 
deflection are compared in figure 10 for two polar angles and two exhaust 
velocities in both the flyover and 85 ° sideline plane. The scatter among 
normalized data points near peak amplitude in the flyover plane was -2.5 dB 
for different polar angles at the same velocity and *4 dB for different vel- 
ocities at the same angle. This scatter greatly exceeds the approximately 
* 3/4 dB day-to-day repeatability of these data. The scatter was caused by 
the previously mentioned spectrum irregularities. All of the noise prediction 
methods use a smooth normalized spectrum or a sum of smooth spectra. There- 
fore the existence of irregular spectra causes an inherent uncertainty in the 
predictions. There was no systematic effect of polar angle or exhaust veloc- 
ity on normalized spectrum. All three predictions generally matched the fly- 
over spectra for Strouhal numbers from 0.2 to 2 and overestimated the levels 
at higher Strouhal numbers. 

A different situation occurred for the takeoff flap setting at 85 ° 
azimuth angle. As shown in the lower part of figure 10, the normalized levels 
at large Strouhal numbers were higher at the larger exhaust velocity. This 
difference occurred because levels for large Strouhal numbers varied with ex- 
haust velocity to at least the eighth power while those near peak amplitude 
followed approximately a 6.5 power variation. At Strouhal numbers larger than 
2, corresponding to highly weighted frequencies for perceived noise at full 
scale, the ANOP prediction matched the data for the higher exhaust velocity 
and rearward direction. The GEIAC prediction was closer to data for the lower 
velocity, and the UTRC prediction generally was between the two. 

Slotless Version of Three -Flap Wing 

Measured directivity in the flyover plane for an UIW slotless wing having 
a lower surface contour tangent to that of the three-flap configuration at 
approach flap deflection is plotted in figure 11. The data are compared with 
predictions by the UTRC and GELAC methods and with the MSA Lewis method of 
references 24 and 25 . This MSA method represents the OASPL noise radiation 
from slotless wings as a sum of two terms. One term is the noise measured 
for a jet impinging against a large flat surface at the same deflection angle 
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and impingement velocity, scaled from the data of reference 3. Because these 
data were influenced by ground reflection, they were decreased 2 dB for 
comparison with these free-field slotless wing data. The other term repre- 
sents trailing edge noise. It uses measured variations of deflected-jet width, 
boundary layer thickness, and maximum velocity at the trailing edge within 
an equation having the correct functional dependence. The GELAC hand calcu- 
lation method of reference 29 predicts that a slotless UTW configuration will 
be 3 dB quieter than the double or triple slotted equivalent shape. 

Additional comparisons of the NASA Lewis and UTRC methods are made in the 
section entitled "Additional Slotless Wings". 

Measured directivity shapes changed from a broad lift -dipole sort of 
pattern at low exhaust velocities, with peak amplitude normal to the deflec- 
ted flap region, to nearly constant amplitude below the wing and flap at 
large velocities. The UTRC method matches these shapes and levels, although 
it generally predicts too low a noise level for directions above the deflec- 
ted aft surface. For directions below the wing and flap, agreement ranges 
from 2 dB underestimate at the lowest exhaust velocity to 2 dB overestimate 
at the highest, levels calculated by the MSA Lewis method of references 24 
and 25 agree with the UTRC predictions near 100° polar angle where both are 
dominated by deflected-jet noise. This MSA Lewis method matches the data 
within 2 to 3 dB over the range of polar angles from 20° to 120° for which 
the method applies. Note that because the test model is relatively small, it 
was necessary to use the flow-field data of reference 24 for the same model 
size rather than that of reference 25 for a geometrically larger model. The 
constant -shape GELAC curve tended to be about 5 dB above the data, and did 
not match measured shapes for the higher velocities. 

Free-field spectra measured at 90° polar angle are compared in figure 
12 with those calculated by the UTRC method. Irregularities, including a 
strong minimum in the I /3 octave band centered at 63 O Hiz, dominate the low- 
frequency portion of the measured spectra. These spectrum irregularities were 
especially strong near 40° polar angle, which is not an obvious angle for noise 
generation or reflection for this slotless wing. Above I 6 OO Hz these spectra 
at 90 ° angle contain mild oscillations of about 2 dB half -amplitude about the 
calculated smooth curves. Oscillations were shown in reference I 9 to be larger 
for this slotless wing than for the three-flap slotted wing at either deflec- 
tion. Possible causes of the nonsmooth spectra were discussed therein. Note 
that according to the data of reference 33 , this slotless wing was too many 
nozzle diameters downstream to produce noise by acoustic feedback between the 
nozzle and deflected solid surface. The measured spectrum oscillations strongly 
affect the statistical accuracy of measixred OASPL. Therefore the normalized 
1/3 octave spectra for this slotless wing are less reliable than those for the 
triple slotted flaps. 
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Calculated and measured normalized l/3 octave spectra for the slotless 
wing are shown in figure 13. The spectrum measured at the lower exhaust 
velocity and 70° polar angle contains the largest scatter. Measured normal- 
ized levels for this spectrum and Strouhal numbers larger than 2 are about 5 
dB below those for the other spectra. All the calculated normalized spectra 
lie 2 tb 3 dB above those other measured spectra for Strouhal numbers larger 
than one. As with all other normalized spectra shown except those for the 
flyover plane and takeoff flap deflection, the AKOP method markedly under- 
predicts measured levels for Strouhal numbers less than 0.2. This may be 
important for predictions of airframe structure acoustic fatigue. 


Comparison With Designated USB Data 


QCSEE USB, Takeoff 

Acoustic and aerodynamic results for I/II.5 and I/28 scale models of the 
QCSEE USB configuration were presented by NASA in reference 22. Acoustic 
data were compared by MSA with predictions made at UTRC before the data were 
available. These comparisons are plotted herein. Calculated and measured 
OASEL directivities in the flyover plane are compared in figure l4 for the 
larger model. Measured directivities clearly changed shape with exhaust ve- 
locity. Amplitudes near 90° polar angle were predicted by the UTRC method 
within 2 dB, and the directivity shapes were closely matched. In contrast, 
the ANOP method matched the measured shapes only at low exhaust velocities 
(not used at takeoff) and was about 4 dB low near 90° polar angle. Ampli- 
tudes near 90° angle, as calculated by the GELAC method of reference 30, were 
no worse than those from the ANOP method, but directivity shapes were in poor 
agreement with data. This GELAC method contains several noise components 
which vary with nozzle velocity raised to the eighth power. For USB config- 
urations they represent noise from the exhaust jet downstream of the trailing 
edge, the wall jet on the wing and flap upper surface, and the canted nozzle 
exhaust impinging against the wing upper surface. However, the sum of these 
calculated terms is 5 to 10 dB below what would be needed to achieve reason- 
able agreement with the data. As mentioned in the preceding section entitled 
"EBF Noise Prediction Methods", even the noise from an isolated exhaust 
nozzle would be calculated as 7 dB less than that given by the NASA ANOP 
method of reference 16. 
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Effects of polar angle on calculated and measured sideline directivity- 
are given in figure 15 for 63 ° and 85 ° sideline angles. Results are shown 
for the nominal flap length and for two other flap lengths to be discussed 
later. Measured reductions be-tween 75° and 105° polar angle were about 6 dB 
and 9 dB for 63 ° and 85 ° sideline angle, respectively. These reductions were 
predicted by the UTRC method and greatly undeip>redicted by the ANOP and GELAC 
Measured reductions between 75° and 105° polar angle were about 6 dB and 9 dB 
for 63 ° and 85 ° sideline angle, respectively. These reductions were predic- 
ted by the UTRC method and greatly underpredicted by the ANOP and GELAC 
methods. The ANOP prediction for USB is defined only for 90 ° polar angle. 
Details of the measured sideline effect in the aft quadrant were poorly pre- 
dicted by the UTRC method. As with the comparison for UTW configurations, 
maximum noise reduction occurred behind 90 ° polar angle but was predicted to 
occur at or ahead of that angle. Increased OASPL at 63 ° sideline angle, at 
a polar angle within the deflected jet exhaust refraction valley in the fly- 
over plane, was larger than predicted and occurred at 20° larger polar angle. 
This good qualitative agreement, but poor agreement in details, probably 
results from the incorrect assumption that quadrupole noise radiation from a 
USB slot nozzle is axisymmetric about the deflected jet centerline. 

Effects of polar angle and exhaust velocity on calculated and measured 
normalized spectra for the I/II.5 scale model are compared in figure 16. Data 
are shown for 60° and 120° polar angles at about I50 and 220 m/sec exhaust 
velocities. Measured normalized spectra for the flyover plane and large 
Strouhal numbers are highest for the aft polar angle at both velocities. 

They are about 8 dB lower for the forward angle and lower velocity. This 
change in spectrum shape and level is predicted by the UTRC method. Although 
the GELAC method gave a poor prediction of OASPL at 120° polar angle, it 
correctly predicted that noise at this direction would be dominated by 
quadrupole noise. The resulting calculated normalized spectra also match the 
high-frequency portion of the data. The ANOP normalized spectrum is inde- 
pendent of polar angle and exhaust velocity; it would generally match the 
data (not shown) for 90° polar angle and both velocities. Measured normal- 
ized spectra at large Strouhal nimbers were greatly affected by exhaust 
velocity at 60° polar angle but not at 120° angle. For 60° angle. Increasing 
the exhaust velocity greatly increases the calculated amount of quadrupole 
noise and therefore the spectrum amplitudes at high Strouhal numbers. It 
causes only moderate increases of surface-radiated noise which dominates the 
peak amplitudes and OASPL. The AWOF method overpredicts measured normalized 
levels by about 3 to 11 dB at high Strouhal numbers. Normalized spectra 
given by the GELAC method were a few dB closer to data than those from the 
UTRC method for Strouhal numbers larger than 5. 
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Three difference flap lengths had been tested with the 1 /II .5 scale 
model at 218 m/sec exhaust velocity. The effect of flap length on directiv- 
ity and spectrum in the flyover plane is shown in figure 17. As can be seen 
from the sketch at the top of this figure, the short flap was only slightly 
shorter than the nominal length but the long flap was considerably longer. 

The short flap was found to cause 3 to 4 dB increase, and the long flap about 
2 dB decrease, of OASPL relative to that for the nominal flap. In contrast, 
the UTRC method predicted no significant change with the short flap and about 
2 dB increase with the long flap. The UTRC method tends to predict an in- 
crease of surface -radiated noise due to increased flap length, and a decrease 
of both quadrupole and surface -radiated noise caused by viscous decay of ex- 
haust velocity at the flap trailing edge. Data given in figure 4 of refer- 
ence 22 showed 10 percent lower peak velocity at the trailing edge of the 
long flap. This velocity difference would be expected to cause 3 <3B noise 
decrease below UTRC calculated levels, resulting in a 2 dB overestimate as 
with the nominal flap. Negligible viscous decay was predicted for the short 
and nominal flap lengths, and very little for the long flap. The UTRC 
method therefore could predict the measured effect of increased flap length 
on OASPL for this configuration if measured velocity at the trailing edge 
could be supplied as input, rather than calculated as one portion of that 
method. , 


Measured normalized I /3 octave spectra for the three flap lengths at 
105 ° polar angle are plotted in the lower part of figure 17. Also shown is 
the normalized spectrum calculated by the UTRC method, which was the same for 
these flaps. This calculation was within 2 dB of data for the nominal and 
long flap. The normalized spectrum measured with the short flap was 2 to 3 
dB higher than the other data near peak amplitude and the same increment 
lower at high frequencies. That is, for frequencies above 435 Hz full scale 
(5000 Hz model scale) which strongly affect annoyance -weighted noise, actual 
measured spectra for the short and nominal flaps agreed within 2 dB with each 
other and within about 3 dB with the UTRC prediction. Large differences be- 
tween 1/3 octave spectra for these two flap lengths were concentrated below 
about 110 Hz full scale ( 125 O Hz model scale), near peak amplitude of the 
spectra. Calculated full-scale anno 3 rance levels therefore would underpredict 
data by about 3 dB for the small and nominal flaps. For the long flap, they 
would be about 5 dB too high. 

All three flap lengths had been tested at 63 ° sideline angle in addition 
to the flyover plane. As shown in figure 15, there was no significant effect 
of flap length on sideline noise reduction. This result is in agreement with 
predictions by all methods. The measured 5 dB reduction at 90 ° polar angle 
was predicted within 1 dB by the UTRC method and underpredicted by the other 
methods . 
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I I I 


Q.CSEE USB. Approach 

Calculated and measured OASEL directivities in the flyover plane for the 
nominal flap length and approach flap deflection are compared in figure l8. 

As with the takeoff configuration, the measured directivity pattern changed 
shape as exhaust velocity was increased. It varied from relatively constant 
amplitude in the forward quadrant at small exhaust velocities to a shape 
which resembled jet noise rotated through the flow deflection angle at large 
exhaust velocities. The ANOP method correctly predicted the measured levels 
ahead of 60° polar angle but increasingly underestimated peak levels at 90° 
as exhaust velocity was increased. The UTRC method correctly predicted the 
change of shape but tended to overestimate by several dB the peak noise for 
directions below the deflected jet. For the typical critical approach case 
of 90° polar angle and 189 m/sec exhaust velocity, measured OASPL was about 
2 dB above the ANOP prediction and an equal amount below the UTRC predic- 
tion. The GELAC prediction for the sum of trailing edge noise and low-ampli- 
tude quadrupole noise terms differed greatly from the data. Underestimates 
of 10 to 13 dB generally occurred near 90° polar angle. 


Effects of polar angle on calculated and measured sideline directivity 
are shown in figure 19 for the approach configuration and 190 m/sec exhaust 
velocity. The measured 5 dB maximum reduction at 90° polar angle and ^3° 
sideline angle was predicted within 1 dB by the UTRC method and underestima- 
ted 2 dB by the ANOP method. Reductions calculated by the GELAC method de- 
crease to zero for directions near the deflected trailing edge, where calcu- 
lated quadrupole noise components are much larger than the calculated trail- 
ing edge noise. Measured maximum noise reduction was only about 1 dB 
larger at 85° than at 63° sideline angle, contrary to the 4 dB increase cal- 
culated by the UTRC method. However, the ANOP method incorrectly predicts 
an increase of OASPL for that change of sideline angle. This relatively poor 
prediction by the UTRC method of sideline noise reduction at approach flap 
deflection and 85° sideline angle may be associated with the underprediction 
of quadrupole noise at polar angles above the deflected flap in the flyover 
plane (figure I8). The slot jet should be loudest at directions normal to the 
narrow side of the jet. For this large sideline angle and large flap deflec- 
tion, the microphone at 90° polar angle is not shielded from line-of -sight 
view of the exhaust nozzle. Any underestimate of exhaust jet noise above the 
wing would cause an overestimate of sideline noise reduction. 


Calculated and measured normalized I/3 octave spectra are compared in 
figure 20 for noise radiated by the QCSEE USB approach configuration. Re- 
sults are shown for 60° and 105° polar angles and nominal 150 and 220 m/sec 
exhaust velocities in the flyover plane, and 60° and 120° polar angles in the 
63° sideline plane. In the flyover plane, the UTRC method predicted a small 
range of variation for normalized spectra. These calculated spectra were 
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close to that given by the AKOP method. In contrast, the GELAC method pre- 
dicted a larger range of spectrum shapes that was in closer agreement with 
data. The same differences occurred between predictions for the 63° sideline 
plane. Spectra measured in this plane were more closely predicted by the 
UTRC method at large and small Strouhal numbers. 

TF-34 Scale Model 

This model, shown in figure 4b, was 1/18.5 the size of a large-scale USB 
model having a 4:1 canted slot nozzle and 40° deflection short flap, tested 
with the mixed exhaust of a TF-34 turbofan engine. Free -field data for this 
geometrically similar scale model had been compared in reference 22 with data 
for that large-scale model, taken from reference 34. OASFL directivity at 
high subsonic exhaust velocities was characterized by a relatively strong 
peak of noise at directions just below the deflected jet. Data from previous 
tests of large- and small-scale USB models had not contained this feature. 

It was not apparent whether this localized apparent quadrupole noise was 
peculiar to the nozzle and flap geometry or was somehow associated with use 
of a real turbofan engine as an air supply. Both this scale model and the 
TF-34 engine installation had uniform exhaust velocities. However, the ex- 
haust of an engine would be hotter and more turbulent than that obtained with 
unheated compressed air, and could have an incompletely mixe'd high-velocity 
core. This small model of the large configuration was tested by NASA to re- 
solve this question. It was shown in reference 22 that the model data, 
scaled to the large configuration, did reproduce the measured directivity and 
spectra. Data were obtained only for the flyover plane. 

Measured and calculated OASPL directivities are compared in figure 28 
for five exhaust velocities. At the lowest velocity (II5 m/sec) the measured 
directivity had no abrupt peak and was closely matched by the ANOP and UTRC 
methods. Increasing the exhaust velocity caused a large increase of noise 
radiation at 100° to 130° polar angles, as with the data of reference 34. 

The AKOP method did not predict this change. The UTRC method predicted the 
qualitative increase of OASPL but underestimated its magnitude. However, 
OASPL at 90° polar angle was predicted within 2 dB by that method for all ex- 
haust velocities. The GELAC method of reference 30 again greatly underesti- 
mated the data for all but the lowest exhaust velocity and less than 90° 
polar angle, as if the quadrupole terms were too low. 

Calculated and measured normalized spectra are compared in figure 22 for 
66°, 90°, and 120° polar angles at II6 and 239 m/sec exhaust velocities. At 
both velocities the UTRC method predicts a narrower range of spectrum shapes 
than the GELAC method and lower normalized amplitudes at high frequencies 
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than the AKOP method. The UTRC method clearly gave the closest prediction 
of data above 10 hHz model frequency (5OO Hz large-scale frequency). This 
would be the dominant region for predicting annoyance -weighted noise levels 
of the large-scale configuration tested with the TF-34 engine. The measured 
rapid decay of normalized amplitudes below 1 kHz model frequency (50 Hz 
large-scale frequency) was not predicted by any of the three methods . The 
AKOP method does predict this type of rapid decay but was 5 to 10 dB above 
the measured levels. 

USB Vane Deflector 


This configuration, shown in figure 4c and described in reference 28, 
has a circular nozzle at moderate height above the wing. A vane deflector, 
mounted on a pivot located above and downstream of the nozzle, forced the 
exhaust jet down against the wing upper surface for powered-lift flight. Ex- 
ternal airflow would pass between the wing upper surface and the exhaust jet. 
During cruise the deflector would be retracted and stowed away, reducing 
aerodynamic friction drag relative to conventional USB (figure 4a). Addition- 
al acoustic data, not given in reference 28, were provided by MSA for use in 
this evaluation of EBF noise prediction methods . 

Calculated and measured OASPL directivities in the flyover plane for 
this configuration at takeoff flap deflection are compared in figure 23 for 
three exhaust velocities. The AKOP method predicts the general level but 
not the shape of these data. The GELAC method matches the measured levels of 
OASPL near 90° polar angle but gives a very poor prediction of directivity 
shape. The UTRC method matches the general shape but is about 7 dB too large 
in amplitude. This large error was caused by the calculated lift fluctuation 
noise component associated with the large chord. However, measured noise did 
not exceed that for the TF-34 scale model which had the same nozzle equiva- 
lent diameter but 0.425 times this chord. Unpublished velocity distributions 
measured for this vane deflector show that the exhaust jet was spread over a 
very large spanvise extent relative to that for USB slot nozzles. Thus the flow 
field achieved with a vane deflector was not typical of that for which the pre- 
diction method is based. 

Calculated and measured sideline directivities for this takeoff config- 
uration are compared in figure 24. At 60° sideline angle (figure 24a), the 
measured reductions relative to those for the same polar angle and the fly- 
over plane were closely predicted by the UTRC method. Maximum measured re- 
ductions of about 5 at two exhaust velocities were about twice the reduc- 
tion predicted by the AKOP method. This close agreement with the noise 
increments predicted by the UTRC method illustrates that the dominant noise 
radiation process below the exhaust jet at both 0° and 6o° sideline angles 
must be surface -radiated noise. 
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In contrast, noise radiation measured at 85 ° sideline angle was stronger 
than that for 60° sideline angle. This result is predicted by the UTRC 
method only when the nozzle exit plane and exhaust -deflecting surfaces can 
be viewed above the wing surface. It is likely that the high -mounted vane 
deflector was not shielded by the limited-span wing model for any polar 
angle directions in the 85 ° sideline plane. Calculated and measured absolute 
values of OASPL at that sideline angle are compared in figure 24b . All three 
methods matched the general level of the data, with the UTRC method being 
about 3 <iB low. It is possible that deflection of the exhaust Jet by an ex- 
ternal circular vane causes more noise than deflection by a nozzle roof in- 
clined at the same angle, as was assumed in the calculation. 

Calculated and measured OASPL directivities in the flyover plane for the 
approach configuration are compared in figure 25 for two velocities . The 
ANOP and GELAC methods generally match these measured levels. The UTRC 
method is about 8 dB too high near 90° polar angle and about 3 dB too high 
at forward positions. That is, the underestimate of surface -radiated noise 
associated with the highly deflected flap was not as severe as the underesti- 
mate of noise associated with the undeflected part of the wing. 

Calculated and measured normalized I /3 octave spectra in the flyover 
plane for the USB vane deflector model are compared in figure 26. Data for 
takeoff flap and vane position, shown in figure 26a, had one broad peak at 
the relatively low Strouhal number of 0.1 and another broad peak at Strouhal 
numbers near 2 . The ANOP method matched the peak at the low Strouhal num- 
ber. However, it and the other two methods predicted a second peak at 
Strouhal numbers from 0.2 to 0.5 where the data had a local minimum. Both 
the ANOP and UTRC methods generally matched the data for Strouhal numbers 
larger than 1. The GELAC method predicted a wide envelope of normalized 
spectra for this range of polar angles. For 120° polar angle the OASPL cal- 
culated by the GELAC method was dominated by noise attributed to direct 
radiation from the wake downstream of the trailing edge. This noise compo- 
nent decays slowly with increasing Strouhal number at this measurement direc- 
tion. Thus the upper dotted line at large Strouhal numbers, which lies 
furthest above the average data, corresponds to the open triangle data sym- 
bols which lie below the average data. 

Measured normalized spectra for the approach configuration, shown in 
figure 26 b, had sharper peaks than those measured for takeoff. The high-fre- 
quency peak, centered at a Strouhal number of 2, extended over about an 
octave of frequency. It protruded more than 6 dB above the remainder of the 
spectrum at the lower exhaust velocity. All of the noise prediction methods 
gave smooth curves which did not reproduce this spectrum peak. The UTRC and 
ANOP methods predicted the general level of data. The GELAC method again 
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predicted a much larger range of normalized spectrum shape than was measured. 
These 80° and 100° polar angles are relatively close to the deflected exhaust 
Jet downstream direction. For Strouhal numbers from 1 to 10, noise calcula- 
ted by the GELAC method for 80° polar angle was predicted to be dominated by 
quadrupole noise from the wake downstream of the trailing edge. At this mod- 
erate angle from the deflected Jet downstream direction, the resulting cal- 
culated spectrum had a large decay rate. For 100° polar angle the calculated 
quadrupole noise radiated from the flap upper surface wall Jet, and refracted 
around the trailing edge, was calculated to dominate at large Strouhal num- 
bers. This normalized spectrum has a small decay rate and produced the 
upper dash curve. For Strouhal numbers near 4, this upper curve is about 8 
dB higher than the 100° polar angle (triangle) data symbols. This strong 
sensitivity of high-frequency spectrum shape to polar angle, at directions 
within about 60° from the deflected Jet downstream direction, is a property 
of predictions by the GELAC USB method of reference 30. The predicted strong 
variation is independent of errors in absolute level of quadrupole noise, and 
is contrary to the data. 

Normalized spectra for the takeoff configuration at 90° polar angle and 
both 60° and 85 ° sideline angle are shown in figure 27. At 60° sideline 
angle, shown in figure 27a, the data for high Strouhal numbers were bracketed 
by the ANOP and UTRC methods. Each method was about 3 dB from the data. For 
85 ° sideline angle, shown in figiure 27b, the ANOP and GELAC methods generally 
were within 3 dB of the data. The UTRC method gave worst agreement with 
data for this configuration and sideline angle. 


Discussion of Error 

When EBF noise data are applied to prediction of full-scale flyover 
noise, calculated levels of perceived noise level (PNL) are dominated by 
the noise radiated at 90° polar angle. A rough estimate of error incurred 
by use of each prediction method was obtained from the difference between pre- 
dicted and measured OASPL for this overhead position. These comparisons 
were examined only at the second highest exhaust velocity for which directiv- 
ity data are presented. Error in predicting annoyance-weighted noise was 
obtained by comparing PNL for scaled predicted and measured spectra at this 
direction and velocity. Model linear dimensions and far-field distance 
were each multiplied by 10 in order to weigh the high-frequency portions of 
the measured spectra in a manner similar to that for PNL at full scale. 
Atmospheric attenuation over the increased path length was included in these 
PNL predictions. The resulting prediction errors for each configuration, 
mean error, and range about that mean for 50^ confidence (0.6T times the 
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standard deviation) are shown in TABLE III for the UTRC, ABOP, and GELAC pre- 
diction methods. These calculations of PNL, and the statistical analysis, 
were conducted hy NASA Lewis Research Center using spectra predicted as part 
of this contract effort. 

TABLE III - OASPL AND PNL PREDICTION ERROR FOR 90° FLYOVER POSITION 

Notes: Errors = Prediction-Data, PNL calculated for 10 times model scale. 

Comparisons for the test exhaust velocity closest to 225 m/sec. 



OASPL Error 

, dB 

PNL 

Error, 

PNdB 

Configuration 

UTRC 

ANOP 

GELAC 

UTRC 

ANOP 

GELAC 

3 -Flap UTW, t/o 

-3.5 

-2.7 

-6.7 

-6.6 

-4.1 

-8.4 

APP 

-1.5 

-2.5 

+0.6 

-0.7 

-1.8 

+0.5 

Slotless 3 -Flap 

+1.5 

NA 

+3.0 

+2.0 

NA 

+2.7 

QCSEE USB, t/O 

+1.5 

-3.2 

-4.9 

+1.1 

-5.2 

-5.7 

APP 

+1.9 

-2.5 

-10.4 

+4.6 

-1.2 

-9.4 

TF-3^ USB 

-1.6 

-h.2 

-5.8 

-3.2 

-6.9 

-6.7 

Vane USB, T/O 

+9.0 

+0.1 

-2.3 

+8.0 

+2.0 

+0.2 

APP 

+7.1 

-0.8 

-0.1 

+11.2 

+2.3 

-3.2 

Mean Error 

+1.8 

-2.3 

-3.3 

+2 .0 

-2.1 

-3.7 

50')(i Confidence 

±2.9 

±1.0 

±3.0 

±3.9 

±2.3 

±3.0 


Both the ANOP and UTRC methods have mean errors of about 2 dB in both 
OASPL and PNL. Mean error of the GELAC method was about 50^ larger. The 
range of scatter about these mean errors (assuming a Gaussian distribution) 
was about half as large for the ANOP method as for the other two methods . 

The UTRC method, best predicted the shapes of the polar and azimuthal 
directivity, but it failed to predict OASPL and PNL at the flyover position 
significantly better than the earlier empirical ANOP method. There are three 
key reasons for this latter result: (l) calculated noise levels are sensi- 

tive to local mean velocity and turbulence level, which are crudely repre- 
sented in the UTRC method, (2) by summing several components whose levels are 
maximum at directions determined by the flap geometry, it is easy to miss the 
level at any one direction, and (3) the ANOP method was based on data corre- 
lations for the specific 90° direction angle at which this comparison was 
made. The resulting ANOP predictions, at this 90° direction angle and 
exhaust pressure ratios near 1.^, apparently are at least as accurate as 
those from the UTRC and GELAC methods which are based on summations of 
separately calculated noise components. If the vane USB configuration had 
not been included in this comparison, the UTRC method would have achieved 
mean errors smaller in magnitude than 0.5 dB, with 50^ confidence levels of 
^1.5 dB OASPL and ^ 2. J PNdB. These predictions clearly were better than 
those of the other methods . 
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One possible alternate method for UIW and USB noise prediction would 
regard measured noise amplitudes and spectra as a sum of two independent 
components: surface -radiated noise and quadrupole noise. Quadrupole noise 

as inferred from the data would be subtracted from those data to find an 
amplitude at 90° direction, directivity shape, and spectrum shape associated . 

with surface-radiated noise. These quantities would vary with configuration 
type (UTW or USB) and flap deflection angle but would be independent of the 
detailed configuration geometry, as with the current ANOP method. Use of 
these two components, each optimized to give best results at the 90° flyover 
direction, should give better prediction of amplitudes, spectra, and general 

directivity shapes for this and other directions. This possible new method 
would require less computation effort than the noise component methods. It 
could be used with UTW mixer nozzle configurations and slotless wing config- 
urations. Use of this empirical approach would give up all pretense of 
describing EBF noise as a sum of simple basic noise processes for which the 
prediction accuracy is limited only by the accuracy of estimated mean and 
fluctuating local flow properties . 

Comparisons for Additional EBF Configurations 
Additional Slotless Wings 

The method presented herein for calculating noise radiated by slotless 
wings had been strongly influenced by data for the UTW slotless version of a 
triple -slotted wing (reference 19)* That paper had included comparisons with 
the calculated sum of a deflected-jet quadrupole noise component and a trail- 
ing edge noise component with amplitude matched with data for the upper for- 
ward quadrant. That sum greatly underpredicted the noise measured below the 
wing at low exhaust velocities. The NASA Lewis method of references 24 and 
25 used essentially these same two components. Attention was confined to the 
limited range of polar angles below the wing and flap. The deflected -jet noise 
data of reference 3 , used in reference 25 for representing quadrupole noise, 
is similarly limited because half the range of polar angles was shielded be- 
hind the large flat surface. Also, data are given therein only for 15°j 30°, 
60°, and 90° deflection without a simple method for interpolating to inter- 
mediate angles. The UTRC method uses an empirical modification to the 
accepted prediction (reference 16) of noise radiated by an isolated jet. 
Therefore it can be readily applied for all deflection angles and exhaust ve- 
locities, and includes noise radiated above the deflected exhaust jet down- 
stream of the slotless flap trailing edge. 

Calculated and measured OASEL directivities in the flyover plane are 
plotted in figure 28 for the large double -slotted wing of reference 25 with 
the slots closed by plug fairings. The calculated curves designated NASA were 
taken from reference 25; they closely match the data. The UTRC calculated 
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curves did not match the fine detail of measured directivity near the deflec- 
ted jet but generally were within 2 dB of the data. Therefore the UTRC 
method for slotless wings gives acceptable prediction of these directivity 
data. 


Calculated and measured directivities in the flyover plane are ccsnpared 
in figure 29 for the UTW slotless wing of reference 26. This configuration 
matched the length and contour of a double-slotted wing at approach deflec- 
tion; it was shorter than the slotless version of the triple -slotted wing. 
Data were presented for two subsonic and (not shown) one supersonic exhaust 
velocity. These data have not been corrected to free field, so all of the 
predicted levels have been increased 2 dB to account approximately for ground 
reflection. 

For the higher exhaust velocity, bhe NASA Lewis method closely matched 
the measured trends and levels at directions below the wing, where this 
method applies. Note that if that method had been assumed to apply at all 
directions, a broad peak of trailing-edge noise centered at 300° polar angle 
would be predicted to occur. Maximum amplitude of that peak would exceed 
the levels calculated by this NASA Lewis method for 20° polar angle. This 
result is inherent in the analytical description and is contrary to the mea- 
sured directivity shapes . The UTRC method predicted the general level but 
not the precise shape of data below the wing. It matched the shape measured 
above the wing, although levels were underpredicted 3 to 4 dB at most direc- 
tions. For the lower exhaust velocity the UTRC method closely matched the 
measured shape and level below the wing and the NASA method was low. The 
GELAC method was about 5 dB high for both velocities. 

Measurements of OASPL directivity and sound power spectra had been 
presented in reference 35 for a UTW slotless wing at several flap lengths. 

The test model had a straight section corresponding to an undeflected wing, 
followed by a circular arc with 45° turning angle. This curved portion ex- 
tended below the nozzle centerline. Straight flap sections with lengths of 
0.75j 3*75j and 12.75 diameters were added to the model. Measured directiv- 
ity in the flyover plane at a jet exhaust Mach number of 0.84 (exhaust veloc- 
ity 270 m/sec) was given in figure 13 of reference 35 for these three flap 
lengths. These data are plotted in figure 30, with the flaps denoted as 
short, medium, and long. 

Calculations by the NASA method of reference 2h required scaling or 
interpolation of the jet impingement noise data of reference 3 to this de- 
flection angle. Those data have a more sharply peaked directivity shape for 
30° than for 60° deflection. The slotless wing data of figure 30 for 45° de- 
flection also have a sharply peaked directivity. Calculations by the NASA 
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method used an interpolation between the noise data of reference 3 for 30° 
and 60° deflection to obtain a prediction for the required 45° deflection. 
Neither the slotless wing data of reference 35 nor the jet impingement noise 
data of reference 3 were corrected to free field. Values for jet exhaust 
width, boundary layer thickness, and maximijm velocity at the trailing edge 
were taken from data of reference 24 scaled for variations with flap length. 

Increased flap length causes increased exhaust cross section area, increased 
boundary layer thickness, and decreased maximum velocity, producing different 
calculated noise for the short and medium flaps. The UTRC method predicted a 
small decrease of noise as flap length was increased from short to medium. 
Noise levels calculated by the UTRC method were increased 2 dB to account for 
the measurements being conducted above a reflecting surface. Calculated 
directivity curves are not plotted for the GELAC method of reference 29. The 
normalized directivity curves given in figure 5 “19 of that report change 
shape drastically between 20° and 60° flap deflection, if the maximum rela- 
tive amplitude shown in that figure for 60° was assumed to apply for 45° de- 
flection, the resulting broad peaks of the calculated curves would match the 
levels of data plotted in figure 30 for the medium flap length at 80° and 
200° polar angles. It would not give the sharp peaks of noise measured near 
the deflected exhaust jet. Because all flap lengths have the same impinge- 
ment-point location, they would all have the same calculated noise levels. 

As shown in the upper part of figure 30, the NASA method matched the 
measured OASPL shape below the wing for the short and mediiom flap. It was 
about 3 dB below data for the short flap and 1 dB above data for the medium 
flap. The UTRC method did not reproduce the measured sharp peak at 110° 
and 120° polar angle, 15° to 25° below the deflected jet, and was 4 to 8 dB 
below data. The UTRC method underestimated angular extent of the noise peak 
measured above the deflected jet but correctly predicted the measured levels 
near 270° polar angle. Levels measured with the short flap generally were 
about 5 dB louder than with the medium flap. Neither method predicted this 
difference . 

Calculated and measured directivities for the long flap are shown in the 
lower part of figure 30. The NASA method predicted the measured sharply 
peaked shape and measured levels . The UTRC method poorly predicted the mea- 
sured directivity shape for polar angles below the wing. A calculated curve 
is not shown for the GELAC method because this geometry is beyond the range of 
variables used in that method. 
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UTW Mixer Nozzle 


Noise radiation data for the nominal half -scale UTW mixer nozzle config- 
uration of reference 37 were compared with those calculated by the method 
given herein. Although the only UTW EBF airplane now flying (McDonnell- 
Douglas YC-15) uses mixer nozzles, no other method has been published for 
predicting the noise of such installations. This method calculates UTW sur- 
face-radiated noise for an exhaust velocity equal to the average of the two 
highest peak velocities in the isolated mixer nozzle’s measured velocity pro- 
file. Measured noise for the isolated mixer nozzle is increased correspond- 
ing to the calculated effect of flow deflection by the flaps at the local 
flow velocity (reference 18). For this configuration, maximum local velocity 
at the impingement point was 0.64 times the jet exhaust velocity. The 
resulting increases were 0.1 dB for takeoff and 0.6 dB for approach deflec- 
tion. This measured, slightly increased noise associated with the mixer 
nozzle flow field was rotated through the flow deflection angle. For polar 
angles above the wing and deflected flaps, it was added directly to the cal- 
culated surface-radiated noise. For polar angles below the wing and deflec- 
ted flap, it was increased 3 dB and added to the calculated surface -radiated 
noise . 

These tests were not conducted in a free -field environment and are 
affected by ground reflections. Acoustic calibration of the test area have 
shown that wave cancellation tended to occur in the range from 400 to 63O Hz 
center frequency. This frequency region, and frequencies from 800 to 2000 Hz 
where ground reflection caused about I.5 dB increase, generally dominated the 
measured OASPL. Tabulated values of OASPL, plotted herein, are believed to 
be about 1.5 dB too large. 

Measured OASPL directivities in the flyover plane, and those calculated 
by the method given herein, are compared in figure 31 for both takeoff and 
approach flap deflections. Results are shown for nominal exhaust pressure 
ratios of 1.2, 1.4, and I.7 at both deflections. Calculated levels and di- 
rectivity shapes closely match the data. At directions near 30° above and 
below the deflected jet, the calculated levels are dominated by measured 
noise from the mixer nozzle as modified in the manner described above. 

Because data for the UTW configuration and the isolated mixer nozzle contain 
the same ground reflection effect, the relative difference between calculated 
curves and measured data symbols is unaffected by ground reflection at these 
angles. Further away from the deflected jet, calculated levels are primarily 
surface -radiated noise. Maximum OASPL was approximately independent of flap 
deflection because it was dominated by noise associated with the mixer nozzle 
rather than the flaps. Increasing the pressure ratio caused the measured, 
nearly flat directivity shape below the wing to develop a peak near the de- 
flected jet. This change was correctly predicted. 
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For direction angles of practical interest, the portion of the frequency 
spectrum which would dominate annoyance -weighted noise would be predicted to 
be associated with the mixer nozzle. Surface -radiated noise would be signi- 
ficant only at lower frequencies. Measured spectra therefore are not com- 
pared with predictions. Such con 5 )arisons were given in reference l8. The 
best evaluation of this noise prediction method for UTW mixer nozzles would 
be comparisons with flyover noise data for the USAF YC-15 Advanced Medium 
STOL Transport. This comparison would have to include measured, rather than 
calculated, forward flight effects on noise from the mizer nozzles. 

Engine in Front of Wing 

Calculated and measured OASPL directivities for this configuration at 
two exhaust velocities for takeoff and approach flap deflection are shown in 
figure 32. Data for this small model (reference 37) were not measvired under 
free-field conditions and are plotted 3 dB below the tabulated levels for 
this comparison. Measured levels were unaffected by axial position of the 
wing leading edge relative to the nozzle exit plane, and this parameter does 
not affect the noise prediction. As compared with UTW and USB configurations, 
the measured directivities are relatively unaffected by flap deflection. 
Calculated noise levels beneath the wing were underestimated 2 to 4 dB at 
the higher velocity but were closely predicted at all directions at the lower 
velocity. 


Calculated and measured spectra at the directions for maximum OASPL are 
compared in figure 33. These spectra are characterized by a 12 dB per octave 
decay at high frequencies. This decay is more rapid than that associated 
with either UTW or USB. Measured spectra generally were closely predicted 
but were underestimated at low frequencies which dominated the contributions 
to OASPL. 


Noise Predictions for Full-Scale QCSEE Configurations 
Zero Forward Speed 

Noise radiation caused by the presence of the wing and trailing edge 
flaps was calculated for the full-scale QCSEE engine, wing, and flap instal- 
lations. Both the UTW and USB configurations were represented at the engine 
exhaust velocities specified as the design takeoff and landing conditions . 
NASA tests of these configurations at zero forward speed are scheduled to be 
run within the next two years. Calculations discussed in this subsection 
were conducted for a 100 m (305 ft) far -field distance and zero forward 
speed. Configuration geometries supplied by NASA are shown in figure 34. 


40 



The UTW configuration had a double -slotted flap with relatively large forward 
flap chord. The aft flap panel extends relatively far into the exhaust Jet 
at takeoff flap deflection. Therefore it is likely that the predictions 
given herein will underestimate the flap-radiated noise for this condition, 
as occurred with the UTW triple -slotted flap model previously discussed. 

Model tests of this QCSEE configuration are recommended, to determine whether 
the increased noise does occur . The high bypass ratio engine had coaxial fan 
and core exhaust Jets. Equivalent exhaust velocity Ve was defined in terms 
of the fan exhaust velocity Vj and exhaust area Ap, and core exhaust velocity 
Vq and exhaust area Aq, in the same manner as with the NASA ANOP method of 
reference 27 . That is, (Ap+Aq)Ve^ = ApVp^ + AqVq^. Other velocity weighting 

functions for defining an equivalent exhaust velocity, such as a mass flow 
weighted velocity (references 22 and 3 ^) , could also have been used. Deter- 
mination of the correct weighting fxmction for two-stream UTW installations 
also requires additional model-scale tests. The sum of fan and core exhaust 
area was I.96 m^ (21.1 ft^) corresponding to I.58I m (5.I88 ft) diameter. 
Velocities specified by NASA for UTW takeoff and approach are: 

Equivalent 

Fan Velocity, Core Velocity, Velocity, 



m/sec 

(ft/sec) 

m/sec 

(ft/sec) 

m/sec 

(ft/sec 

Takeoff 

204 

(670) 

245 

(803) 

215 

( 704 ) 

Approach 

146 

( 480 ) 

194 

( 636 ) 

161 

(528) 


The USB configuration had an internal mixer nozzle. Fully mixed exhaust 
velocities were specified as 220 ra/sec (722 ft/sec) for takeoff and I90 m/sec 
(623 ft /sec) for approach. Nozzle equivalent diameter was 1.491 m ( 4.893 
ft). Atmospheric properties were taken as those for standard sea level. 

Because actual exhaust velocities for these planned full scale tests may 
differ from those specified here, and other definitions of equivalent veloci- 
ty may be examined, it is necessary to predict the effect of small changes in 
exhaust velocity. This was done by regarding OASPL at each polar angle, 
sideline angle, and flap deflection as varying with equivalent exhaust veloci- 
ty raised to some exponent n. OASPL was calculated to two decimal places for 
the nominal and 1.0233 times nominal velocity. Ten times the difference 
between the two values of OASPL was then equal to the velocity exponent. 

The calculated OASPL directivities in the flyover plane and at 63° and 
85° sideline angles are plotted in figure 35 for UTW takeoff. Also shown are 
the calculated variations of velocity exponent with polar angle for these 
three sideline angles. The velocity exponent at 90° polar angle is largest 
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for 85 ° sideline angle because calculated noise at this direction is dominated 
by quadrupole noise. Calculated OASPL directivities at these sideline angles 
for UTW approach, and the associated velocity exponents, are plotted in figure 
36 . For both takeoff and approach, the calculated OASPL in the flyover plane 
increases by 1 to 3 dB as polar angle is increased from upstream toward the 
deflected jet. Predicted velocity exponents have maxim'um values of about 9 
near the edge of the jet exhaust refraction region. These values are some- 
what larger than tend to be actually measured. Calculated 1/3 octave spectra 
at 90 ° polar angle and the three sideline angles are plotted in figure 37 for 
both takeoff and approach. Maximum 1/3 octave SPL is predicted to occur below 
50 Hz center frequency, the lowest frequency used in many methods for predic- 
tion of annoyance-weighted noise, for five of these six spectra. If OASPL is 
measured as the simi of 1/3 octave SPL's for center frequencies from 50 to 
10,000 Hz, this quantity would be about 3 dB less than the actual OASPL plotted 
in figures 35 and 36 . 

Calculated OASPL directivities and velocity exponents for the QCSEE USB 
configuration at take off and approach conditions are plotted in figures 38 
and 39- OASPL in the flyover plane is predicted to have very little varia- 
tion with polar angle at directions beneath the deflected exhaust jet. Again, 
the calculated velocity exponents near the edge of the exhaust jet refraction 
region seem unrealistically high. Calculated levels of OASPL at 90° polar 
angle in the flyover plane are between 102 and 103 dB for both UTW and USB at 
takeoff. They are about 97 and 98 dB for the two configurations at approach. 
The specified geometries and exhaust velocities therefore are well balanced 
to provide calculated conditions of about equal noise for both UTW and USB. 
Calculated 1/3 octave spectra at 90° polar angle for USB takeoff and approach 
are plotted in figure ^tO. As with the UTW spectra, only about half the OASPL 
was predicted to be radiated above 50 Hz center frequency. 

Effects of Forward Flight 

Calculations were conducted of the effects of forward flight on spectra 
at a position 100 m (328 ft) directly beneath the QCSEE configurations. The 
flight speed was taken as 4l m/sec (80 knots), the QCSEE nominal design con- 
dition for both takeoff and approach. As was previously mentioned, this cal- 
culation method predicts spectra that would be measured in a coordinate 
system which is fixed. relative to the airframe. Such predictions are 
appropriate for comparison with data for models tested in open jets or 
acoustic wind tunnels. They include the effects of a reduction in noise 
source strength caused by changes in turbulence level of the exhaust jet. 

They do not include the changes in directivity pattern caused by motion of 
the noise sources relative to the atmosphere. This change would divide the 
mean square acoustic pressure by the quantity 1-Mj.cos 0 raised to some posi- 
tive exponent, where Mp is the flight Mach number. However, EBF noise 



annoyance is predicted to be largest for polar angles 0 near 90° > and EBF 
flight Mach numbers for takeoff and landing are near 0.2. Thus the omitted 
correction is unimportant for conditions of practical importance. Predic- 
tions for a coordinate system fixed relative to the aircraft can be converted 

to those for ground-fixed microphones by dividing mean square acoustic pres- 
sure by 1-I^cos G and Doppler-shifting the frequency. This correction also is 
small for cases of practical interest. Calculated spectra are shown for the 
flyover position, at which the omitted factors do not change the predicted 
noise. 

The calculated effect of 4l m/sec flight velocity on flyover spectra 
for the QCSEE UTW takeoff and approach configurations is shown in figure 4l. 
Each predicted spectrum is decreased by nearly a constant increment at all 
frequencies. The decrease is slightly larger for approach, with its larger 
ratio of flight velocity to exhaust velocity, than for takeoff. In contrast, 
the same comparison for the QCSEE USB configuration (figure 42) shows only 

about 1 dB predicted noise reduction at frequencies of interest for noise 
annoyance. This small effect occurs because forward flight is predicted not 
only to reduce the spectrum amplitudes but to shift each amplitude to a 
higher frequency. It was shown in references l8 and 20 that these qualita- 
tive trends have been reported for tests of EBF configurations in acoustic 
wind tunnels. 


CONCLUSIONS 


1. The UTRC method best predicted the variations of EBF noise amplitude 
with polar and azimuthal angle. The UTRC, ANOP, and GELAC methods 
adequately predicted normalized spectrum shapes. 

2. All three methods poorly predicted OASPL at the flyover position for 
some configurations. The UTRC method failed badly for a USB circular 
nozzle with vane deflector at both takeoff and approach flap deflections, 
but generally was closest to data for the more conventional UTW and 

USB Installations. This discrepancy probably was caused by differences 
between the actual and calculated or assumed local flow-field proper- 
ties. The consequence of this poor agreement for two of the eight 
example cases was that the ANOP method gave more acciirate predictions of 
average PNL at 90° direction than did the UTRC or GELAC methods. 
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BECOMMEM)ATIONS FOR FUTURE RESEARCH 


Flow -field turbulence properties such as turbulence streamwise and 
transverse integral scale length and turbulence convective velocity, along 
with mean and rms fluctuating velocity, should be measured for simple slotted 

UIW and unslotted USB configurations. These data should be utilized along 
with available theories for calculating the spectrum of lift force fluctua- 
tion and noise radiation. For USB, it may be necessary to develop a numeri- 
cal simulation for the pressure field induced by a random distribution of 
discrete vortices representing the wall -jet boundary layer and upper shear 
layer . 


Noise measurements should be obtained for USB configurations having 
conventional nozzles and at least a factor of 2 variation in the ratio of 
upper surface flow length to nozzle diameter. These data would be utilized 
to determine whether discrepancies between measured noise radiation for USB 
models and levels predicted by the method given herein were caused by incor- 
rect prediction of the effect of this parameter. If so, the prediction 
method should be modified. 

Noise measurements should be obtained with a scale model of the QCSEE 
UTW configuration at takeoff flap deflection. Current EBF noise prediction 
techniques poorly predict data for this type of flap position relative to 
the exhaust nozzle . 

Flyover noise data should be obtained and compared with predictions by 
this and other methods. Ideally, these tests should be obtained with a 
powered sailplane or some other type of aircraft having highly suppressed 
engine noise. Tests could be conducted with the two USAF Advanced Medium 
STOL Transport configurations (a USB and a mixer nozzle UTW) if it is pre- 
dicted that flight conditions exist for which EBF noise exceeds noise 
radiated directly from the propulsive systems. 
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The following three annual reports and two final reports were prepared 
under this Contract and published as NASA Contractor Reports. 

Fink, M. R. : Investigation of Scrubbing and Impingement Noise. 

NASA CR-134762, Feb. 1975- 

Fink, M. R. : Prediction of Externally Blown Flap Noise and 

Turbomachinery Strut Noise. NASA CR-134883, Aug. 1975* 

Fink, M. R.: Additional Studies of Externally Blown Flap Noise. 

NASA CR-135096, Aug. 1976. 

Fink, M. R.: A Method for Calculating Strut and Splitter Plate Noise 

in Exit Ducts - Theory and Verification. NASA CR-2955, 1978. 

Fink, M. R.: A Method for Calculating Externally Blown Flap Noise. 

NASA CR-2954, 1978. 

The following AIAA papers, and subsequent publications of some of these 
papers, provided wider distribution of major results obtained under this 
Contract . 

Fink, M. R.: Experimental Evaluation of Theories for Trailing Edge and 

Incidence Fluctuation Noise. AIAA J., Vol. 13, No. 11, Nov. 1975, 
pp 1472-1477. Also, Paper 75-206, AIAA, Jan. 1975. 

Fink, M. R. : Scrubbing Noise of Externally Blown Flaps. Progress in 

Astronautics and Aeronautics, Vol. 45, Aeroacoustics : STOL Noise, 

Airframe and Airfoil Noise, M.I.T. Press, Cambridge, Ma., pp 3 "25, 1976 
Also, Paper 75-469, AIAA, Mar. 1975- 

Fink, M. R. and Olsen, W. A.: Comparison of Predictions and Under-the- 

¥ing EBF Noise Data. Paper 76-50I, AIAA, July 1976. 

Fink, M. R.: Forward Plight Effects on EBP Noise. Paper 77"13l4, AIAA 

Oct. 1977. 



APPENDIX A : LIST OF SYMBOLS 


a Speed of sound, m/sec 
c Total wing and flap chord, m 
Cji Chord of nth flap panel, m 

D Nozzle exit diameter or hydraulic diameter, m 
f One -third octave center frequency, Hz 

h Average distance from flap panel to assumed vortex trajectory, m 
Kn Amplitude function for fluctuating lift noise of nth flap panel 
Me Convective Mach number for Jet, 0.62 Mj 

Mj Jet exhaust Mach number relative to ambient speed of sound 

Mjj Jet exhaust Mach number relative to Jet speed of sound 

Pref Reference acoustic pressure, 2 x 10“5 N/m^ 
r Far -fie Id distance, m 

St Strouhal number, fD/V 

Ui Maximum exhaust velocity at impingement distance, m/sec 
V Nozzle exhaust velocity, m/sec 

X Axial distance from nozzle exit to impingement point, m 

Deflection angle of last flap segment, deg 

0 Polar angle relative to nozzle upstream direction, deg 

9n Polar angle relative to upstream direction along chord of nth flap 

segment, deg 

p Density, kg/m^ 

0 Azimuth angle relative to flyover plane , deg 

Subscripts 

a Ambient atmosphere 

ISA International standard atmosphere 

J Jet exhaust 

L Lift fluctuation noise 

n Nth flap segment 
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The following AIAA, paper, prepared by NASA Lewis Research Center, 
compared NASA data with predictions calculated under this Contract. 

Olsen, W. A., Burns, R., and Groesbeck, D.: Flap Noise and Aerodynamic 

Results for Model QCSEE Over-the-Wing Configurations. Paper 77 "23, 

AIAA, Jan. 1977. 

The following AIAA paper and publication, while, not conducted under this 
Contract, described a direct extension of contract results to an additional 
practical application. 

Fink, M. R.: Approximate Prediction of Airframe Noise. J. Aircraft, 

Vol. 13, No. 11, Nov. 1976, pp 833-834. Paper 76-526, AIAA, July I976. 
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APEENDIX C : COMPUTER PROGRAM FOR CALCULATING EBF NOISE 


General Description 

This digital computer program, written in FORTEIAN IV, predicts 
externally blown flap noise that would be measured in the free field at 
points on a sphere centered at the nozzle exit. These points are equally 
spaced in polar angle, in planes of designated sideline angle. The reference 
angle direction (zero polar angle for all azimuth angles) is forward along 
the nozzle centerline. If atmospheric properties are not specified as input 
but allowed to remain at their sea -level standard default values, all linear 
dimensions should have the dimensions of meters and velocities should be in- 
put as meters per second. Comment statements are placed throughout the pro- 
gram listing to describe the purpose of each portion of the program and to 
define the program variables. The program is relatively small, requiring 
less than 12K of computer memory. 

Input quantities include geometric properties of the EBF configuration 
as sketched in figure ^3. For under -the -wing (UTW',CONFIG=l) and engine-in- 
front-of-the-wing (lPW,CuNFIG=3) installations, these properties include the 
coordinates of the wing and flap leading edges in a coordinate system center- 
ed at the nozzle exit. Here, positive X is downstream and positive Y is up- 
ward, normal to the wing plane. Other geometric properties are the wing in- 
cidence relative to the nozzle centerline, deflection of each flap segment 
relative to the wing, number of flap segments, chord of the last flap segment, 
nozzle exit diameter, and far -field radius. Geometric lengths required for 
upper surface blowing (USB,C0NFIG=2) installations are the coordinates of the 
wing leading edge and flap trailing edge, nozzle exit hydraulic diameter, and 
far-field radius. Geometric angles for these installations are the nozzle 
roof angle (also called cant angle or kickdown angle) relative to the nozzle 
upstream centerline, wing incidence relative to the nozzle upstream center - 
line, and flap deflectic» relative to the wing. An USB configuration without 
a flap is described as having zero flap deflection. 

Two special cases are under -the -wing configurations having slotless wings 
or mixer nozzles. Slotless wings are represented by C0NFIG=1 and NFLAP=0 in- 
put. Wing geometry is input as the leading edge coordinates XW, YW, the de- 
flected trailing edge coordinates X(2), Y(2), incidence DELW of the forward 
undeflected portion of the slotless wing relative to the nozzle centerline, 
and deflection DEL(i) of the aft portion relative to the forward portion. 
Slotted under -the -wing configurations having mixer nozzles are designated 
separately (C0NFIG=4) . They require the same input as C0NFIG=1 plus the iso- 
lated-nozzle centerline velocity ratio at the impingement distance, input as 
ROOF. 
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A flow chart for the calculation process is shown as figure 44. The 
first portion of the computer program calculates various geometric quantities 
such as chord lengths. For UTW and USB, it also calculates distances from 
the nozzle exit plane to the impingement point and the trailing edge, average 
distance from each wing and flap segment to the assumed far edge of the Jet 
shear layer, and resulting spanwise -coherent lift force fluctuation. Axial 
distances downstream of the nozzle exit plane are utilized to calculate the 
ratio of local maximum velocity to nozzle axial velocity. The program vari- 
ables XW and XL(n) are the ratios of wing and Nth flap segment maximum fluc- 
tuating-lift mean square acoustic pressure to reference pressure squared. 

(This noise component has also been called scrubbing noise or inflow noise.) 
The program variable XTE is the ratio of maximum trailing edge noise mean 
square acoustic pressure to reference pressure squared, calculated only for 
the most rearward trailing edge. Calculated quadrupole noise from an iso- 
lated exhaust Jet, in the direction perpendicular to the Jet centerline, is 
adjusted for local velocity ratio and Jet deflection angle to yield the 
acoustic pressure ratios for quadrupole noise of a Jet deflected by UTW flaps 
or an USB nozzle and wing assembly, and of such a Jet downstream of the 
trailing edge. These pressure ratios apply for the direction perpendicular 
to the deflected Jet. 

Next, the noise components are calculated at equally spaced increments 
DELTH of polar angle at up to ten designated sideline angles PHI. For each 
polar angle, the components of fluctuating lift noise from the wing and flaps 
are computed and summed. Similarly, the viewing angle relative to the center - 
line of the deflected Jet must be computed so that quadrupole noise can be 
determined. Calculations of quadrupole noise pressure ratio GJET must take 
into account whether the combination of polar angle and sideline angle yields 
a direction within the Jet refraction region and, for USB, whether the direc- 
tion is shielded from noise generated near the nozzle exit. Individual l/3 
octave spectra for the different noise components are computed and added to 
obtain the spectrum and OASPL of the UTW or USB configuration. 

In contrast, quadrupole noise from engine-in-front -of -the-wing 
installations is taken as that from an undeflected isolated exhaust Jet. Lift 
fluctuation noise is calculated for loading distributions which become acous- 
tically noncompact along the scrubbed span. The l/3 octave spectra of noise 
from the wing and flap is calculated taking this into account, and this spec- 
trum is summed to obtain overall noise. 
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Input Variables 


The following is a list and definition of the input variables 
default values. 

and their 
Default 

Eroeram Symbol 

Definition 

Value 

CONFIG 

Integer equal to 1 for UTW, 2 for USB, 3 for IFW, 
or 4 for UTW with mixer nozzle 

1 

NFLAP 

Integer, number of flaps 

1 

NPHI 

Integer, number of sideline angles (maximum 
allowable number is 10) 

1 

V 

Exhaust Velocity- 

100. 

VINF 

Flight Velocity- 

0. 

D 

Nozzle hydraulic diameter 

1. 

R 

Far -field radius 

1. 

THL 

Lower value for polar angle, deg 

0. 

THU 

Upper value for polar angle, deg 

H 

CO 

o 

• 

DELTH 

Increment between successive polar angles, deg 

10. 

HH(1) 

First polar angle, deg 

0. 

FHI(2) 

• 

Second polar angle, deg 

0. 

• 

• 

PHI (10) 

Tenth polar angle, deg 

0. 

FL 

Lower limit for l/3 octave center frequency, Hz 

50, 

FU 

Upper limit for l/3 octave center frequency, Hz 

• 1E5 

CA 

Ambient speed of sound 

340.3 

RHOA 

Ambient air density 

1.225 
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Default 


Program Symbol Definition Value 

PREF Reference acoustic pressure .2E-6 

DELW Wing deflection relative to nozzle, deg 0. 

DEL(i) Deflection of first flap relative to wing, deg 0. 

DEL(2) Deflection of second flap relative to wing, deg 45. 

• 

DEL(KFLA.P) Deflection of last flap relative to nozzle, deg 45. 

EOOF USB nozzle roof angle, deg, and mixer nozzle 0. 

velocity ratio at impingement distance 

XW Axial position of wing leading edge 0. 

X(l) Axial position of first flap leading edge 1. 

• 

X(NFLAP) Axial position of last flap leading edge 1. 

YW Vertical position of wing leading edge 0. 

Y(i) Vertical position of first flap leading edge 0. 

Y(NFLAP) Vertical position of last flap leading edge 0. 

CLAST Chord of last flap 1. 


The input variables CONFIG, NFLAP, and NPHE must be integers. This program 
can be utilized with English-system units (ft, ft/sec) for lengths and veloc- 
ities if CA, RHOA, and PREF are supplied in English units. For standard 
atmosphere, these are set equal to lll6., .2377E-2, and .41773E-6, respective- 
ly. Data input uses the standard NAMELIST format. A title card, prepared 
for each run, provides a printed title to identify each specific configura- 
tion. 
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Test Case 


The following cards provide a test case for this computer program. All 
of these cards are punched starting in column 2. The first card provides a 
title to identify the print -out. Input dimensions are in meters and exhaust 
velocities are in meters per second, so the default values for standard sea- 
level air density and reference acoustic pressure are used. This test had 
been performed during cold weather, and speed of sound is input for the 
actual air temperature. The inputs for direction angle will cause calcula- 
tions to be performed for polar angles from 60° to 120° in 20° increments, at 
0° and 63° sideline angles, for both exhaust velocities. 

QCSEE USB 1/11.5 SCALE MODEL, TAKEOFF FLAP SETTING 
$INPUT 

C0NFIG=2, NPHI=2,PHI(1)=0.,PHI(2)=63. 

THL=60 . , THU=120 . , DELTH=20 . , FL=50 . , FU=20000 . 

DELW=5 . ,XW=- .20 , IN=- .024 , D= . 123 ,.R=6 . 1 
DEL(1)=30.,X(2)=.42,Y(2)=-.122,CA=330. 

V=191. 

$END 

$INPUT 

V=219. 

^EX^ 

$INPUT 

IEND=2 

$END 
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->41 


1 * 

2 * 

3* 

4* 

5* 

6 * 

7* 

8 * 

9* 

10 * 

li* 

12 * 

13* 

14* 

l5* 

16 * 

17* 

18« 

19* 

20 * 

n: 

23* 
24* 
25* 
26* 
27* 
28* 
29* 
30* 
31* 
32* 
33* 
34* 
35* 
36* 
37* 
38* 
39* 
4 0* 
41 * 
42* 
4 3* 
44* 
45* 
46* 
47* 
48* 
49* 


c*** this program calcul'ates far-fielo externally blown flap noise for 

C*** UNOER-THE-WING (C0NFI6=1), UPPER SURFACE BLOWING ICONFIG = 2) AND 
C*** IN-FRONT-OF-THE-WIN6 (CONFIG = 3> ENGINE INSTALLATIONS* OASPL AND SPECTRA 
C*** ARE CALCULATED FOR THE FREE FIELD AT POINTS ON A SPHERE CENTERED AT THE 
C*** EXHAUST nozzle AND MOVING WITH THE AIRCRAFT AT A LOW SUBSONIC FLIGHT 
C*** SPEED. DEFAULT ATMOSPHERIC PROPERTIES AND REFERENCE ACOUSTIC PRESSURE ARE 
C*** IN METRIC SYSTEM UNITS CORRESPONDING TO INPUT DISTANCES IN METER^S AND 
C*** EXHAUST VELOCITY AND FLIGHT VELOCITY IN METERS PER SECOND. 

c*** mixer nozzle 

REAL KW,K,M,L .LTE.MJ 
DIMENSION DTHETA(37),DEL(4) 

DIMENSION FB{4D),X(4),Y(4),PHI{1D>,CC(41,HL(4J,HT(4),K(4»,XL(3), 

1 SPL(4D,37), 0ASPLC37) 

DATA V, VINF,D,R,THL,THU,DELTH,PHI(1) ,PHI(2) ,PHI(3>,PHI(4),PHI(5», 

1 PHI(6)fPHI (7) ,PHI(8>,PHI(9), PHKID), CA , RHOA , PREF , DEL W , 

2 DEL(1»,DEL(2) ,DEL(3),0EL(4),R00F,FL,FU,XW,X(1>,X(2),X(3), 

3 Y(1),Y(2),Y(3», CLAST, YW/ 

4 IOD.,0.,1. ,1. .0.,180. ,1D.,1D*D. , 340.3,1.225, .00002, C., 

5 D.,3*45.,0»,5u*,1DDDD*,D*,3*1*,3*0*,1.,0./ 
namelist/input/ IEND, CONFIG, NFLAP.NPHI, 

1 V, VINF,D,R,THL,THU,DELTH,PHI (1) ,PHI(2) ,PHI(3) ,PHI(4),PHI(5), 

1 PHI (6),PHI (7) ,PHI(8),PHII9>, PHKIO), CA , RHOA , PREF , DEL W , 

2 DEL (1),DEL (2> ,0EL(3),DEL(4 ) , ROOF , FL ,FU, XW, X ( U , X< 2 ) , X ( 3 ) , 

3 YW,Y(1),Y(2),Y(3), CLAST 

DIMENSION TITLE(13) 

READ{5,85DI (TITLE(I),I=1,13) 

FB(1> = ID. 

FB(2J = 12.5 

FB(3) - 16. 

PB(4> = 20. 

FBI5) =25. 

FB(6) : 31.5 
FB(7) = 40. 

FB(8) = 50. 

FB(9) = 63. 

FB( 10)= 80. 

DO 4 1=1,10 
FB(I+10) = 10.* FB(I) 

FB(I+20) = 100.* FB (I ) 

4 FB(I+30) = lOOC.* FB(I) 

6 REA0(5, INPUT) 

NTHETA= ( THU-THL )/DELTH ♦1.00001 
IFCNTHETA.GT. 37. 1 ) N THETA =37 . 0000 1 
URITE(6,7) 

7 FORMAT(lHl) 

IFdEND.NE .0) GO TO 2000 
WRITE(6, INPUT) 

IF=1 


Program Listing 



\J1 

CO 


5n* 


DO 10 1=1,38 

r>i* 


IF( FL.LT .F B ( I ) ) 00 TO 20 



IF = 1 

5 3> 

10 

CONTINUL 

54* 

20 

LF=38 

55* 


DO 30 I=i,33 

56* 


J = 39-1 

57* 


IF( FU.6T.FB ( J ) ) GO TO 4Q 

58* 


LF=J 

59* 

30 

CONTINUE 

6 0* 

40 

IF(ir.LT.LF) GO TO 50 

61* 


LF=IF 

62* 

50 

CONTINUE 

63* 


OE&^^AD = .01745329 

64* 


IF< XW.LT.O, >XW=0. 

65* 


IF(C0NFIG.EQ.4) GO TO 60 

66* 


IFCCONFIG .NE. 1 ) GO TO 200 

67* 

60 

CONTINUE 

68* 

C ** * 

GE0i3ETRY CALCULATIONS FOR UNDER -T HE -W ING BLOWN FLAPS 

69* 


IF( NFLAP.NE .0 ) GO TO 100 

7 0* 

c*** 

SLOTLESS .UNDE R-T HE- WING CONFIGURATION 

71* 


GO TO 205 

72* 

65 

LTE=X( 1)+CC(1 ) 

73* 


L = LTE 

74* 


C=: 3.*CW 

75* 


IF( Y(2).G£.Q. ) C=CW 

76* 


ARG=4DELW+DEL(1>)*0EGRA0 

77* 


IF( Y(2).LE,0. )L=X(1 )+Y{l)*COS (ARG»/SIN(AR6) 

78* 


IF(Y(l).LE.O.)L-XW*YW/TAN(OELLl 

79* 


GO TO 120 

SO* 

c*** 

SLOTTED FLAPS 

81* 

. 100 

CW = SCRT ( (X ( 1 )-XW ) **2 . + { Y ( 1 )-YW )**2 . ) 

82* 


ARG = DELW*DEGRAD ♦ OEL ( NFL AP )*DEGRAD 

83* 


X(NFLAP+1)= X(NFLAP) + CLAST * CO^IARGJ 

84* 


Y(NFLAP+l)r Y(NFLAP) - CLAST ♦ SlN(ARG) 

85* 


SUM=CW 

86* 


DO 110 N=1»NFLAP 

87* 


CC(N) = SORT( (X(N+1 )-X (N ) )**2 .♦( Y (N + 'l ) -Y (N) )**2.) 

88* 

110 

SUM = SUH+CC(NJ 

89* 


C = SUM 

90* 

c*** 

SHEAR LAYER TRAJECTORY AND FLAP LIFT FLUCTUATIONS 

91* 

120 

HLU = .5*0 +YW 

92* 


HTW = HLW -CW*SIN (DELW*DEGR ad ) 

93* 


IF(HTW.LT,.5*D) HTW =.5*D 

94* 


KW = l.E-7 *C *CW *(HLW**{-2.) + HTW**(-2.) ) 

95* 


IF ( NFLAP.EQ.O ) GO TO 500 

96* 


LTE = X(NFLAP) + CLAST ♦ COS(ARG) 

97* 


L = X(NFLAP) ♦YINFLAP )/T ANt ARG) 

98* 


IF(<Y(NFLAP»-CLAST*SIN(ARG)).GE.O.) L=LTE 

99* 


DO 130 N =1,NFLAP 

100* 


HL(N) = .5*0 +Y(N) 



1 


VJI 

VO 


101 

102 

103 

104 

105 

106 

107 

108 
109 


12,7 

128 

129 

130 

131 

132 

133 

134 

135 


136 


137 

138 

139 

140 

141 

142 


143 

144 

145 

146 
,147 

148 

149 

150 


HT(N) :: HL(N)-CC(N) ♦ SIN ( DELW+DEGRAD +DE L ( N ) *DEGRA 0 ) 

IF(HL f N) ,L T . . 5*0 ) HL(N) ;,5*0 
IF(HT( N) .LT,. 5>»-D ) HT(N) =.5*D 
130 K(NJ r 1 ,E-7*C*CC(N » ♦ ( HL I N ) ♦ ♦ ( -2 . > +HT (N ) ♦♦ ( -2 . ») 

60 TO 500 

200 IF( CONFIG .NE. 2) 60 TO 300 

geometry calculations for upper surface blown T'LAPS 

NFLAPri 

ARG = (DELW+DEL (1) )*,0EGRA0 

205’ CONTINUE “ “ 

OELL-DELW+DEGRAD 

IF(DEL 11 ) .LT. .01 ) X(1)=X(2) 

IF(0EL(1 ) .LT. .01 ) GO TO 210 

X{ U -X (21 -( YW-Y( 2)- (X (2 )-XW > ♦TAN (DELL ) ) / ( TA N ( DELL + DEL ( 1 ) ♦DE6R ADI -T 
lANCDELLM 

210 Y(1)=YW~(X(1)-XW )♦TAN{OELL) 

CW = SORT! (X ( ll-XW )^^2. ♦ ( Y ( 1 ) -YW )* + 2 . I 
CC( 1) - SORT I (X (2)-X( 1 ) y♦♦2. + ( Y (21 -Y ( 1) )♦♦2.), 

IF(C0NFI6 .EQ. 1 IGO TO 65 
C = CW +CC ( 1 J 
L=X(11+CC(1) 

LTE:::L 

KWr2.E-7^C^X( 11/ (D*D) 

K(l)=l .5E-7^C^CC(1)/(D*0) 

GO TO 500 
300 CONTINUE 

C*++ ENGINE IN FRONT OF THE WING 
CROrCLAST 

IF ( NFLAP .EO.O ) GO TO 320 

CW = SQRT( (X( 1 )-XW) ♦♦Z.* (Y (1 l-YW l^+Z. ) 

X(NFLAP+1) = X(NFLAP) + CLAST ♦ COS ( OELW^ DE GRAD+DE L ( NF L AP 1 ♦DEGRA D ) 
Y(NFLAP+1> - Y(NFLAP)' - CLAST ♦ S IN ( DELW + DE GRAD + DEL ( NFL AP ) ♦DEGRA D ) 
DO 310 N=l, NFLAP 

310 CC(N) : SQRT( ( X(N-*-l )-X (N n^^2 . + ( Y ( N + 1 ) -Y ( N I I ♦♦ 2 . ) 

320 IF ( NFLAP .EQ .0 I IUPPER-1 
IF (NFLAP. EQ.l ) IUPPER=2 
IF ( NFLAP .GT .1 ) lUPPER^NFLAP 
M r V/CA 

XJET^ (l.-VINF/V )♦♦6.♦1,E14 ♦ 7 . 5/ ( R + R^ ( 1 . + , 01 ♦K + + 4 , 5 ) ) 

Br .62+M /((I. + ( ,62*M l+^S. I^^O.Z I 
MJ r{ .62^M)+*0.1 
PI =3.14159 

PLK= 1.45E-4 *D *(Pl ♦RH0A*VA+3.* 0/ ( 2 .♦CA^R+PREF ) ) +^2 , 

60 TO 530 
500 M=V/CA 

C*** VELOCITY RATIOS AT IMPINGEMENT POINT AND AT TRAILING EDGE 
IF(C0NFIG.EQ.4) 60 TO 501 

VR=(l.+{(.14+L/n)/SQRT(l.+M)) ++4 . )♦♦ ( -.25) 

VRTE = ( 1. + ( ( . 14*L TE/D)/SORT ( 1 . +M))«^4 . )♦♦ (-.25) 

GO TO 502 



o 

o 


151*- 

152^ 

153^ 

501 

154^ 

155* 

,156^ 

157^ 

502 

158^ 

159^ 

160^ 

161* 

162+ 


163^ 

164^ 

c*** 

16 5 ♦ 

505 

166^ 

506 

167^ 

C*A ♦ . 

168^ 

169^ 

170+ 

c*** 

171 + 

nz* 

173 + 
174^ 
175^ 
lib* 
111* 
178^ 
179 + 

507 

180^ 

c*** 

181^ 
182* 
183 + 
184^ 
185^ 

c*** 

186^ 

510 

187^ 

188^ 

189^ 

190^ 

191^ 

192* 

c*** 

193+ 

194^ 

520 

195^ 

196^ 

197* 

53Q 

198^ 

1994 

ZOO* 



V-V^RCOF 
VRrl . 

VRTE-1. 

CON TI NUE /V j^^2,*c*D*tVR*RH0A*(V*VR CA*PREF+-R ) )**2 . 

IF tCO NF I6.E0.4 JP =2 .♦P 

PPHI - PHKIPHI) ♦OEGRAD 

MAXIMUM IRaTLING EDGE NOISE 

XTE- 1 .E -S*CA=«'P + D*I VR TE/ VR )**'3 . /(C*V ) 

IFt CONF I6.E0. 1 . A NCi.MFLAP.NE .0 »P=. 

XL W - P ♦K W 

IF tNFLAP .EQ .0) GO TO 506 

MAXIMUM LIFT FLUCTUATION NOISE FOR EACH FLAP 
DO 5Q5 N-l.NFLAP 

XLINI=P*i »N1 

NOISE^PERPENDICULAR TO JET WITH VELOCITY EVALUATED AT NOZZLE EXlTf IMPACT 
POINT, OR TRAILING EDGE 
IFJCONFIb.NE.'t) 60 TO 5Q7 

M=.l+H 

CONTINUE 

XJET-n.-VINF/V)**2.*l,Elt|*D*D*M**7. 5/ IR*R* (l. + .01*M**4.5)) 

XM V R — H ♦ VR 

XIH JE T=^ ( 1. - VINF/ V)**2.*l .El 4=*0*D*XMV R**7 .5/(R*R*(l.+ . 01*XHVR ♦♦4. 5 ) 

XMVRTE - M*VRTE 

XTEJET=( I.-VINF/ V)**2.*l ,E1 4*D*D>AXMVRTE + *7, 5/ ( R*R* 1 1 , ♦ . 01 *XM VRTE** 

4.5 n 

IF(C0NF16.NE.2160 TO 510 
INCREASED OUADRUPOLE NOISE 
WING SURFACE 

AMPDEL-1 .+6.*SIN tR00F*0E6RAD )**2. 

XiMJETzX JET*A>4P0EL 
XTEJET=XTEJET*AMPDEL 
6(5 TO 52D 
CONTINUE 

INCREASED OUADRUPOLE NOISE FROM DEFLECTION OF jET BY UTW FLAPS 
AA = U .-VINF/VJ^*4, 

XT- JET^XTEJE T* (1 . + COS ( PPHI » **2 . !♦ (AA+6.*VR**8 .♦SIN (AR6»**2. I /2. 

XIMJET = XIMJET^ tl .♦C0StPPHIJ^^2.}»tAA+6.^SIN(AR61*»2.» /2. 

IF! NFLAP.EQ.OI X IH JET = . 25+ X IMJET 
IFINFLAP.E0.HI X TE JET=. lO^XTEJET 
B = .62^M/lf l. + t .6 2+HI#^5. I )♦♦G.2 
HJ=t .62+Hl^*0.1 

continj E 

DO 999 IPHl^lfNPHI 

WRITE(6,36D> CTITLECIII ),III=1,13 ) 

WRITEI6,B10) PHIIIPHU 
DO 998 ITHETA-1 .NTHETA 

THETA = (THL ♦ (I THE TA-1 ) ♦DEL TH > ♦DEGRAD 


FROM DEFLECTION OF JET BY USB NOZZLE ROOF AND 


CT\ 

H 


201* 

202* 

203* 

204* 

205* 

C *** 

206* 

207* 

208* 

209* 


210* 

211* 

bOO 

212* 

213* 

214* 

215* 

605 

216* 

217* 

218* 


219* 

606 

220* 

607 

221* 

222* 

C*** 

223* 

224* 

225* 

22b* 

227* 

228* 

229* 

230* 

~ 610 

231* 

620 

232* 

233* 

c***- 

234* 

630 

235* 

236* 

640 

237* 

238* 

c*** 

239* 

650 

240* 
24 1* 

660 

242* 

243* 

670 

244* 

245* 

680 

246* 

690 

247* 

248* 

249* 

250* 

i 


DTHtTA(irHETA)-THETA/DEGf^AD 
TRIG z 5IN(THETA)X«SIN (PPHI ) 

ANGF z (COS I A TAN ( rRIG/SCJRT ( 1 . -TRIG * TRIG))))«*2, 

IF(CONFIG .EO. 3) GO TO 740 

CALCULATING SUN OK LIFT FLUCTUATION NOISE COMPONENTS 
IF ( NFL AP .EQ.U ) GO TO 605 

t9T5LHA’fl*9S<OELW«DEbRAD)*SlN( THET A ) C 0 S IP PH I ) +S IN ( DE LW *OEG R AD ) *00 S 

11 IHL \ n ) t 2 m 
SUM zu, 

DO bOD N=1,NFLAP 

SUM-SUM-*-XL(N)*(COS((DELW + DEL(N))*DEGRAO)*SIN(THETA) ♦COST PPHI ) +SI N ( 
1 (DE L W + DEL(N))>^DE GRAD)* COS (THE TAT )*«2. 

GL-Gw^'SUM 


GO TO 607 
GLzl. 

IF( ( THLTA + ARG ) .LE. 3,1416) GO TO 606 
IF( THE TA.GE.6.283) GO TO 606 
GO TO 607 

GL-XLW*(C0S(AR6)*SIN(THETA)*C0S{PPHI)+SIN(ARG)*C0S,( THETA) )**Z, 
C ON T I N UE 

GE z aNGF *XTE *C0S ( (THETA+ARG)*.5 )**2. 

VIEWING ANGLE RELATIVE TO DEFLECTED JET 
ANG JzsiN (PPHI )*SIN(THETA) 

ANGJ =ATAN( ANGJ/ SQR T ( 1 , -ANGJ*ANGJ ) ) 

ANGJ=ATAN (COS ( AN GJ ) **2 . *S IN ( ARG ) / COS ( ARG ) ) 

IF ( ( ( THE T A + ANG J ) *M J ) , LE.2 . 6 1 8 ) GO TO 680 

IF( ( THETA+ANG J ) ,L E. 3. 1 416 ) GO TO 610 

1F( ( ( 6 .28 32-THET A-ANGJ )*MJ J . GE .2 .6 18 ) GO TO 660 

GO TO 69 U 


IF( CONF IG ,EQ . ? ) GO TO 640 

GJETz(xrMJET+XTEJET)*tl,+B*C0S(THETA+ANGJ))**(-3.) 
IF(CONFIG .EQ. 2) GO TO 650 
REFRACTION IN QUAORUPOLE NOISE 

GJETz!5JET*10.**(-1.8*(1.-ABS(6. -(THETA +ANGJ)*1.91*MJ))) 
GO TO 700 


GJETZXTEJET*{ l. + 8*C0S (THETA + ANGJ) )*♦( -3. ) 
IF(C0NFIG.E0.2) GO TO 630 
REFRACTION IN OUADRUPOLE NOISE 


G JETzGJET + 10. ** (- 1 . 8* ( 1 . -A8S ( 6 , - ( 6 . 28 32-THETA 
GO TO 700 

IF(C0NFIG.E0.2) GO TO 620 
GO TO 690 


-ANGJ) *1 .91*HJ ) ) ) 


GJETZXTEJET*( l.+B*COS (THETA+ANGJ) )**(-3.) 

GO TO 700 

IF(C0NFIG.EQ.2) GO TO 67q 

GJET=( XIMJET+XTEJET ) * ( 1 . +B*CO S ( T HE T A+ ANG J ) )*♦( -3.) 

IF(C0NFIG.EQ,2) GO TO 700 

IF ((( THETA+ANGJ ) *HJ ) .LT.2. 618 )G0 TO 70(j 

IF( ( (6.2832-THETA-ANGJ)*MJ).LT.2.618)60 TO 700 

GO TO 650 


251* 
252* 
253* 
254* 
255* 
256 * 
257* 
256* 
259* 
260* 
261* 
262* 
263* 
264* 
265* 
266* 
267* 
268* 
269* 
27 0* 
271* 
2 7 2* 

27 3* 
274* 
275* 
276* 
277* 
278* 
279* 

28 J* 
281* 
232* 
283* 
284* 
285* 
286* 
287* 
288* 
28 9* 

290* 

291* 

292* 

293* 

294* 

295* 

296* 

297* 

298* 

299* 

300* 


700 6JET-GJET*(1.*ISIM«ANGJ)/SIN{AR&))**2.) 

C*** SUM OF NOISE COMPONENTS 
6- 6L + GE + GJET 
OASPL ( ITHETA ) - 10,* AL0G10(6) 

IF ICCNFIG .EQ. 21 GE-GE + 6L 
PLL^Qj, 

00 73Cf I=IF,LF 

C*** CALCULATION OF 1/3 OCTAVE SPECTRA 
S -FB(I)*D/V 

PLJET :: GJET *S**4,/ (8. * (.125 + S** 1 . 3 33 ) ♦♦ 4 . ) 

IF (CaNFIG .EQ . 1 ) GO TO 71Q 
S=S/( 1 .♦ VINF/ V ) 

GO TO 720 

710 PLL = GL *.037 * S**4./ (0.008 +S**Z. 667 1**2. 

720 PLE - GE *.C29 * S**4./ (0.05 +S**1.5 1**4. 

PL = PLL + PLE ♦ PLJET 

730 SPL(I, ITHETA) - 10 . »ALOG 1 0 ( PL ) 

GO TO 998 

C*** ENGINE IN FRONT OF WING 
740 CONTINUE 

C*** QUADRUPOLE NOISE OF UNOEFLECTED JET 
GJET- XJET *( 1 .*6*00$ (THETA ))+*( -3 . ) 

IF( (THETA*MJ) .LE.2.618) GO TO 744 

IF ( THE TA .LE . 3 . 1416) 60 TO 742 

IF( ( ( 6 .2832-THET A )*MJ ) . LE.2 .6 13 ) GO TO 744 

GJET - GJET *10. **(-.9 * ( 1 . -A 3S C 6 . - ( 6 . 23 32-THE TA ) *1 . 9 1 ♦ MJ M ) 

GO TO 744 , 

742 GJET r GJET *lC.**(-.9 *(1.-A8S(6.- THETA*1.91* MJ) ) > 

744 CONTINUE 

C*** LIFT FLUCTUATION NOISE FROM WING AND ALL BUT LAST FLAP, BY METHOD OF 
C*** FILOTAS 
SUMirO . 

DO 300 I-1F,LF 
SUM-0. 

CRD- CW 

DELLr OELW*DEGRAD 
DO 790 jrl,IUPPER 
IF(J.EO.l) GO TO 75C 
CRD :: CC(J-l) 

DELL = (DELW + U E L ( J-1 ) ) * DE GR A D 
750 SC - FB(I) * CRD/V 
SIGMA = PI*SC 
BETA - 2.*CRD/D 
IF ( SIGMA .GE .BETA ) GO TO 770 
AAA = BETA 

760 XLS r (LOGd,? + ( P I * A A A ) * *2 . ) ) / ( L 0 G ( 1 . 2 ) *3 , * ( P I *A A A ) *♦ 2 . ) 

GO TO 780 
770 AAAz SIGMA 
GO TO 76U 

7 80 PLFr SC**3.*(4./ (1. + (.5:j*SI0MA*D/CPD)**2,))*XLS/( 1. + (SIGMA*M)**2.) 


(i J * 

L' .■;>> 
7 

3 [! f) * 
(,!£'=.' 
3 C 7 
3 lit-* 

Z [■ L- 

7 ic.>;- 

31 

3 1 Px- 
3 1 3x> 
3 1 V 
3 3 Sx 
3 1 6 X. 
3 ] 7 >:-• 

3 1 ':■■ * 
2 ;’ux= 
3 2 IX' 
5 2 2 x< 
3 2 3X' 
3 2 *' x< 

3 2 X“ 
326*' 
327X' 
3 26 X' 
3 2 '>x- 
3 3 b*' 
3 5 IX' 
3 3 2 ■■’ 
3 3 3*' 
3 3 X' 
3 3 b X' 
3 3 6 X' 


I L 2 : j. . L r / { V, ' 

;,Ur- rbUr:+ ^'LF ^41. { C';- b ( DLLL ) =-=r. H. (T ht T ^ ) X. Lf:S ( I- Ptll ) ■•• b II- ( DELL ) X-CCS ( IHT T 

1 » . ) ) X'X'/ . 
vc.u C ctM r I r-.'UL 
P L ft - j I ) f i 

5.14-1 OF M.-15-L CC"IOM. MS 
SHM1- S14’ 1 -^FLV 

fLJET= SOFT *■ ( F''- ( I ) 4L / V ) X'X'H . / ( t . X' ( . K ( F L ( I ) X- u / V ) >•■ x - 1 . 3 .5 3 ) ^ * L . ) 

oLii SFH I , 1 I hf T/ ) = ir.,v /'LOClCUf L., -^FLJ^T) 

Ot SF'L ( n hL I f. ) - jL.<= i'LCC-.ir + L-DL I ) 

‘/S'e. crMiF'Ll 

C X- Xt X< 

C CLTFl'I 

c >jt # * >;< 

IFG = L 

k;IjO it = 1 + ]7v ii-r 

IT :: FT Kf (M F L T M n + T 6 ) 

MJTE(e,tM:! < UU LT fi ( I J , I = n , L I ) 

V. F-^ i T £ { 6 , t 2 b ) 

L: 0 K; 1 L i - i F , L F 

10 IL F4,n £ ( t, t. :.L ) F £• ( 1 ) , (SPL.( I , J ■) , jrjl , 1 . 1 ) 

I' R n L ( 6 , L C: ) ( C / £ F l < J ) , J - 1 1 ,L.l ) 

if'C- - IF 6 + 1 

T r ( ( L T -M ||£ T A ) . L T . 2 ) *. F- T It ( C , 7 ) 

IF ( (L T-r,THE TA ! .L T .L ) GO TO IDOC 
9x79 OOMIF'uF, 

0 0 1 0 ' t 
20 rc CONTI f:OL 

blL FCl'MiAT ( // , ?U> , ’ Sif.E Lir.t A OOL E ,F'H T r ’ , F 7 . 2 , 1 X , * OF 0 R [ f S ' it ) 

P21.I fORF’A T( // ' Ml OF*, / ’ CTh Fk£(:’, 1 2 X , • A 1 I MU T H ANGLE, DEG’,/ 

1 • ( H L F- T 2 ) ’ , 1 X . 1 7 E . 0 ) 

82b FOPMAK/) 

l 30 rORUAT ( 2X,F 9 .2, 17F 7.2 ) 

;'.4Lj F0RMAT(/,4X, ’ 0 V f. R A L L ’ , 1 7 F 7 . 2 ) 
hbU F OPMA T ( 1 3A6 ) 
t (. 0 F (. K Fi A T ( U-l 1 , 1 0 X , 


1 3 A f- ) 



Test Case Output 


SINPUT 

lEND 

- 

♦0 

CONFIG 


+ 2 

NFLAP 


+ 1 

NPHI 

z 

+ 2 

V 


.219CQC00E+03 

VINE 


.DOQOQGDDE+DO 

D 

z 

. 1230aOOOE + DO 

R 

z 

,6lCOaOOGE+01 

THL 

z 

.60DG0DG0E+02 

THU 

z 

. 12Q0000CE + 03 

DELT H 

z 

.2D000DDDE+D2 

PHI (1) 

z 

, QUOQQOUCE + aO 

PHI (2) 

z 

,6300UOUnE-*-02 

PHI (3) 

z 

. aOQOOOOOE + OO 

PHI ( 4 ) 

z 

. uCOOODQOE + 00 

PHI (5) 

z 

.aoooQcaoE+QO 

PHI (6) 

z 

. OJOOOOOOE+ DC 

PHI (7) 

z 

,JOuOOOJD£+GO 

PHI (8) 

z 

. OQOODCODE-*-00 

PHI (9) 

z 

. UODOQOJOE +'JD 

PHI (10) 

z 

. OODOOCOOE + 00 

CA 

z 

. 33C000a0E+ 33 

RHOA 

z 

. 1225u00nE + Ql 

PREF 

z 

.2a00G0'JDE-04 

DEL1» 

z 

.50Q0ODC0E+J1 

DEL (1) 

z 

. 3O0QOCDDE + O2 

DEL (2) 

z 

,45b0u00CE+02 

DEL (3) 

z 

. 45UOOOOOE +02 

DEL (4) 

z 

.4500U000E+D2 

ROOF 

z 

. OOCOOGCOE+OQ 

FL 

z 

.5D00UG00E+02 

FU 

z 

•.20C00QOGE + D5 

Xkl 

z 

, 03COGOOnE + CO 

X ( 1 ) 

z 

. 32002793E+C0 

X (2 ) 

— 

.420GUC00E+00 

X ( 3 ) 

z 

, lDOOOOOnE + 01 

YW 

z 

-.2400DCG0E-01 

Y ( 1 ) 

z 

-.5 199881 lE-01 

Y (2 ) 

- 

122CUG0DE + 0D 

Y (3 ) 

z 

. OOOOOOUGE + 00 

CLAST 

z 

. lOuOuDuDE + 01 


SEND 


CCSEE USb l/il.5 SCALE M 0 D EL , T A K E 0 F F FLAP SETTING 
SIDELINE ANGLE, PHI- ,00 DEGREES 


1/3 OE 

CTR FREQ: AZIMUTH ANGLE, DEG 


(HERTZ) 

6 0. 

8 C . 

ICO. 

120, 

SO. 00 

76,40 

75.93 

74.86 

73.04 

6 3 , uO 

79.8 1 

79,34 

78.29 

76.49 

8 0.00 

83.11 

62.65 

81.60 

79 .8 4 

100.00 

35.92 

85.47 

84.45 

8 2.72 

125.00 

88.42 

87.98 

86.98 

85.32 

160.00 

90.74 

9 0.32 

89.37 

87.81 

2L0.00 

92. 39 

9 1.99 

91.10 

89.66 
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25D.00 93.58 93.22 

315.00 94.32 94.02 

400.00 94.60 94.36 

500.00 94.45 94.31 

630.00 93.96 93.93 

800.00 93.15 93.26 

1000.00 92.20 92.44 

1250.00 91.09 91.48 

1600.00 89.74 90.29 

2000.00 88.44 89.12 

2500.00 87.09 87.90 

3150.00 85.65 86.58 

4000.00 84.13 85.19 

5000.00 82.71 83.86 

6300.00 81.22 82.47 

8000.00 79.69 81.03 

10000. CO 78.26 79.68 

12500.00 76.84 78.33 

16000.00 75.28 76.83 

20000.00 73.88 75.48 


92.41 91.13 

93.31 92.25 

93.81 93.01 

93.93 93.40 

93.75 93.51 

93.31 93.36 

92.70 92.98 

91.94 92.43 

90.95 91.63 

89.95 90.77 

88.87 89.81 

87.69 88.73 

86.41 87.54 

85.18 86.37 

83.88 85.13 

82.51 83.82 

81.23 82.57 

79.93 81.31 

78.49 79.90 

77.18 78.62 


OVERALL 104.01 103.98 103.81 103.57 


QCSEE USS 1/11.5 SCALE MODEL , TAKEOFF FLAP SETTING 
sideline ANGLEfPHl= 63.00 DEGREES 


1/3 OB 
CTR FREQ: 


AZIMUTH ANGLE, DEG 

(HERTZ ) 

60. 

80. 

100. 

120. 

50.00 

71.65 

69.52 

68.04 

66.83 

63 . 00 

75.06 

72.93 

71.47 

70.31 

80.00 

78.36 

76.24 

74 . 80 

73.70 

100.00 

8 1.18 

79 . 07 

77.66 

76.65 

125.00 

8 3.68 

8 1.59 

80.22 

79.34 

160.00 

86.01 

83.95 

82.64 

81.95 

200.00 

87.66 

85.64 

84.42 

83.96 

250.00 

88.86 

86.89 

85.79 

85.61 

3 15 .00 

89.62 

87.72 

86.78 

86.96 

4 CO. 00 

89.9 2 

88 . 12 

87.39 

87.97 

SCO. 00 

89.80 

88.12 

87.62 

88.59 

630.00 

89.34 

87.81 

87.57 

88.92 

8 to . CO 

88.58 

87.22 

87.25 

88.95 

i J'wC . 00 

87.67 

86.48 

86. 76 

88.73 

1 2lD. 00 

86.62 

85.59 

86.09 

88.29 

16^0.00 

85.32 

84.47 

85.20 

87.59 

Z 0 uC . 0 0 

84.07 

83.37 

84.27 

86.80 

25uC. 00 

82.76 

82.21 

83.25 

85.89 

3 1 bC. CO 

81.37 

80.95 

82.12 

84.85 

4 0^0. 0-J 

79.90 

79.60 

80.89 

83.70 

5 0 L, 0 . u 0 

78.51 

78.31 

79.69 

82.57 

■ 6360.00 

77.06 

76.96 

78.42 

81.35 

& Li u 0 . 0 0 

75.57 

75.55 

77.08 

80.06 

1 D - j c G . no 

74.17 

74.23 

75.82 

78.83 

125^0.00 

72.77 

72.90 

74.54 

77.58 

16 0 uO .00 

71.24 

71.43 

73.12 

76.18 

20060.00 

69.86 

70.10 

71.83 

74.91 

0 V er all 

99.40 

97.86 

97.62 

98.98 
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SINPUT 



TEND 


♦ 0 

CONF IG 


+ 2 

NFLAP 

z 

+ 0 

NPHI 

3 

+ 2 

V 

3 

.19100000E+03 

VINF 

Z 

.OOQOOOQOE+OO 

D 

z 

• 123D0000E + 00 

R 

3 

.61000D0DE+01 

THL 

3 

.60DOaOCOE+02 

THU 

3 

. 12C00000E + 03 

DELTH 

z 

.2000DOD0E+Q2 

PHI (1) 

z 

.OOOOOOOOE+OO 

PHI (2) 

z 

.6300D0D0E+02 

PHI (3) 

3 

.ooooooooE+oa 

PHI (4) 

3 

.ODDOODDOE+DO 

PHI (5) 

3 

.QOODOOOOE+DO 

PHI (6) 

3 

.□OOCOOOOE+DC 

PHI (7) 

3 

.GOOOODUaE+OO 

PHI (8) 

3 

.ooooaocoE+oc 

PHI (9) 

Z 

.ODOODOODE+OO 

PHI (10 

> = 

.□ODOUOOQE+00 

CA 

3 

. 33ODOD00E^03 

RHOA 

3 

. 1225UD00E + Q1 

PREF 

3 

.20D00000E-04 

DELU 

3 

.50000CQOE+Q1 

DEL (1) 

3 

.3aaDDDUQE+D2 

DEL (2) 

3 

.45D0D0QDE+D2 

DEL (3) 

3 

.4500DQ00E+02 

DEL (4) 

Z 

.450D0000E+02 

ROOF 

z 

.OOOOOOOOE+OO 

FL 

3 

.5000UC00E+02 

FU 

3 

.20000000E+05 

XU 

Z 

-.20QOOOOOE + 00 

X (1 ) 

3 

. lODOOQODE + Dl 

X (2 ) 

z 

.42CguD00E+0D 

X (3 ) 

3 

. lOOOOOOOE+01 

YU 

z 

-.24000000E-01 

Y (1 ) 

3 

.OOOOOOOOE+OO 

Y (2 ) 

3 

-.i22DOOOOE+DO 

Y (3 ) 

Z 

.OOOOOOOOE+OO 

CLAST 

z 

. lOOOOOOOE+01 

SEND 




QCSEE USB 1/11.5 SCALE MO DEL , T A K EO F F FLAP 
SIDELINE ANGLE, PHI= .DO DEGREES 

1/3 OB 


CTR FREQ : 


AZIMUTH ANGLE, DEG 

(HERTZ) 

60. 

80. 

100. 

120. 

50.00 

75.12 

74.55 

73.35 

71.31 

63.00 

78.41 

77.84 

76.65 

74.63 

80.00 

81.55 

80.98 

79.81 

77.81 

100.00 

84.17 

83.61 

8.2.45 

80.49 

125.00 

86.44 

85.89 

84.76 

82.85 

160.00 

88.47 

87.94 

86.84 

85.02 

200.00 

89.83 

89.32 

88.27 

86.55 


SETTING 


66 



250.00 

90.72 

90.24 

89.26 

87 

.68 

315.00 

91.15 

90.72 

89.83 

88 

.45 

4 00.00 

91.14 

90.77 

90.01 

88 

.86 

5 00.00 

90.75 

90.45 

89.84 

88 

.94 

630.00 
8 00. 00 

90.03 

89.83 

89.39 

88 

.75 

89.04 

88.95 

88.70 

88 

• 32 

1000.00 

87.93 

87.96 

87.88 

87 

.74 

1250.00 

86.69 

86.84 

86.94 

87 

.00 

1600.00 

85.22 

85.50 

85.78 

86 

.03 

2000.00 

83.82 

84.22 

84.66 

85 

.06 

2500.00 

82.38 

82.90 

83.47 

84 

.00 

3150.00 

80.87 

81.49 

82.20 

82 

.83 

4000.00 

79.29 

80.02 

80.84 

81 

.58 

5000.00 

77.80 

78.63 

79.56 

80 

.37 

63CO.OO 

76.27 

77.19 

78.21 

79 

.09 

8000. 00 

74.69 

75.70 

76.80 

77 

.74 

10000.00 

73.22 

74.31 

75.48 

76 

.47 

12500.00 

71.77 

72.93 

74.15 

75 

.19 

16000.00 

70. 17 

71.40 

72.69 

73 

.76 

20000.00 

68.75 

70.03 

71.37 

72 

.47 

overall 

100.54 

100.29 

99.74 

98 

.95 


QCSEE USB 1 / 11.5 SCALE HO DEL » TAKEOFF FLAP SETTING 
SIDELINE ANGLE, PHI= 63.00 DEGREES 


1/3 OB 
CTR FREQ: 


AZIMUTH ANGLE, DEG 

(HERTZ) 

60. 

80. 

100. 

120. 

50.00 

70.43 

68. 14 

66.54 

65.00 

63.00 

73.72 

71^44 

69.85 

68.35 

80.00 

76.86 

74.58 

73.01 

71.58 

100.00 

79.49 

77.22 

75.68 

74,33 

125.00 

81.76 

79.51 

78.01 

76.77 

160.00 

83.80 

81.57 

80.13 

79,08 

200.00 

85.16 

82.98 

81.62 

80,78 

250.00 

86.06 

83.93 

82.68 

82.11 

315.00 

86.52 

84.45 

83.35 

83.12 

400.00 

86.53 

84.56 

83.65 

83.80 

5 00.00 

86.17 

84.31 

83.61 

84.13 

630.00 

85.49 

83.78 

83.31 

84. 19 

8 00.00 

84.54 

82.99 

82.77 

83.98 

1000.00 

83.48 

82.08 

82.09 

83.56 

1250.00 

82.30 

81.06 

81.27 

82,96 

1600.00 

80.88 

79.81 

80.23 

82. 12 

2000.00 

79.54 

78.62 

79.20 

81.23 

2500. 00 
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APPENDIX D: NOISE SOURCE LOCATION FOR UPPER SURFACE BLOWING 


Purpose of Experiment 

Several analyses have been presented for the acoustic process by which 
noise is generated by USB externally blown flaps. The physical situation 
along the jet centerline plane is sketched in figure 45* An engine exhaust 
nozzle is located on a wing upper surface. The presence of the wing causes 
observers beneath the wing to be shielded from much of the aft-radiated en- 
gine noise. When the trailing edge flaps are deflected, their upper surfaces 
form a smooth low-curvature surface. The exhaust jet remains attached to 
this curved upper surface and is deflected downward, generating significant 
increases of wing lift at low flight speeds. Broadband noise with a velocity 
exponent less than that for isolated exhaust jet noise is also generated. 

At small distances downstream of the exhaust nozzle, the exhaust jet can 
be regarded as containing a potential core having uniform mean flow proper- 
ties and low turbulence level. A boundary layer exists between the potential 
core and the wing upper suface, and a shear layer exists above the potential 
core. The upper shear layer has a faster growth rate and higher turbulence 
level than the boundary layer. Depending on the configuration, these two 
layers of strong viscous effects may merge upstream of the flap trailing 
edge. The boundary layer can still be denoted as a distinct region because 
of its larger gradients of streamwise mean velocity and smaller normal tur- 
bulence levels, caused by the boundary constraint at a solid surface. Down- 
stream of the trailing edg^e, the upper surface boundary layer becomes a shear 
layer. The large gradient of mean velocity within the boundary layer pro- 
duces a large growth rate and high turbulence levels in the initial part of 
this shear layer. Thus the largest turbulence levels in any portion of the 
exhaust jet near the trailing edge are those of the shear layer downstream of 
the trailing edge and below the exhaust jet. 

One analysis of USB noise, developed by Tam and Reddy (reference 2l), 
assumed that one of the dominant components was quadrupole noise generated 
in the lower shear layer. The calculated levels depend upon measured proper- 
ties of the turbulence and the mean velocity profile. Some USB noise reduc- 
tion concepts tested by Hayden (reference 38 ) have assumed that the noise is 
caused by turbulent eddies convected past the trailing edge and spatially 
very near that edge. Thus the noise is implicitly assumed to depend on tur- 
bulence properties of the wing upper surface boundary layer immediately up- 
stream of the trailing edge. A noise prediction method developed by Fink 
(reference 17) represented USB noise at low exhaust velocities as a sum of 
two simple surface-radiation noise processes: trailing-edge noise and 
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whole-body lift fluctuation noise caused by large-scale structure of the 
exhaust jet as it is convected along the airfoil and past the trailing edge. 
The upper shear layer ahead of the trailing edge, with a combination of large 
turbulence scale length and large turbulence intensity, was regarded as the 
dominant region. It seemed likely that crosscorrelation between local flow 
velocities and far field acoustic pressure should indicate which region, if 
any, contains the dominant noise source. 


Apparatus and Procedure 


Test Apparatus 

Tests were conducted in the anechoic chamber of the UTRC acoustic wind 
tunnel. A 15 cm (6 in.) inside diameter duct brought a regulated muffled 
supply of high-pressure air into the chamber. This air supply duct was con- 
nected to a 4.9 cm (1.925 in.) exit diameter axisymmetric nozzle. A 23 cm 
(9 in.) chord MCA OOI8 airfoil was mounted such that for zero deflection 
relative to the nozzle centerline, the airfoil leading edge was vertical. Its 
closest position to the nozzle was 2.5 cm (l.O in.) downstream and 2.0 cm 
(0.8 in.) to one side of the nozzle lip. The airfoil was rotated about its 
30 percent chord line to 9° deflection, trailing edge toward the exhaust jet. 
This nozzle and airfoil had previously been used in tests reported in refer- 
ences 17j 18, and 20. When positioned in this manner, the aft 70 percent of 
the airfoil surface adjacent to the jet was nearly in line with the nozzle 
lip. A conventional 0.635 cm (l/4 in.) microphone was located 2.13 m (7 ft) 
to the side of the nozzle exit centerline shielded by the- wing, to measure 
far -field acoustic pressure. Thus the test configuration was an USB instal- 
lation rotated 90°. 

When this test program was planned, it was intended to measure the 
exhaust jet flow properties with a miniature hot film gage used in tests 
described in reference 20. However, repeated irregular output spectra 
occurred. A conventional single hot wire gage, which generally is less rugged 
than hot film gages, was therefore used. This hot wire gage was mounted on 
a traverse mechanism which could be manually positioned axially in a plane 
normal to the airfoil chord and containing the nozzle centerline. The hot 
wire could be remotely actuated in horizontal position within this plane. 

The amplified output signal could be connected to circuits which measured 
mean output and rms fluctuation, corresponding to mean and rms axial veloc- 
ity. 
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Procedure 


Par -field acoustic spectra were measured with and without the hot wire 
gage and its support strut extended across the exhaust jet at 100 m/sec (328 
ft/sec) exhaust velocity. The presence of the support strut caused up to 2 
dB noise increase at 2500 Hz frequency and up to 7 dB increase at 6300 Hz. 
This larger frequency was the expected Strouhal frequency for vortex shedding 
by the 0.32 cm (l/8 in.) diameter hot-wire support strut. Output of the far- 
field microphone and the hot wire was therefore sent through matched filters 
adjusted to pass only the frequency region between 100 and 2000 Hz. The re- 
sulting filtered far -field acoustic signal was essentially unaffected by hot- 
wire position. The dominant broadband peak amplitude occurred at about 400 
Hz frequency and was unaffected by the filters. 

The hot wire was traversed across the exhaust jet to obtain profiles of 
mean velocity and rms velocity fluctuation in the direction parallel to the 
nozzle centerline. These traverses were conducted at axial locations half a 
nozzle diameter upstream and downstream of the trailing edge. Positions of 
maximum rms turbulence in the wing boundary layer, lower shear layer, and 
upper shear layer at both axial stations were located during these traverses. 
The hot wire was then positioned at each of these four locations. Far -field 
acoustic spectrum was measured in the frequency range from 100 to 2000 Hz to 
assure that the spectrum was unaffected by the presence of the hot-wire 
support strut. The measured difference between noise signals with and with- 
out the probe was largest at the support probe Strouhal frequency near 10,000 
Hz but was greater than 0.2 dB down to 3150 Hz. Autocorrelations of the re- 
sulting filtered broadband hot-wire and acoustic-pressure signals, and a 
crosscorrelation of those signals,, were obtained. Tests were also conducted 
with the hot wire located at the same vertical distances relative to the air- 
foil surface at axial positions one diameter upstream and downstream of the 
trailing edge. All data were obtained at 100 m/sec exhaust velocity. 

Hot-wire autocorrelation traces were integrated with respect to time to 
obtain the Euler ian integral time scale. This time scale was multiplied by 
the mean velocity at the measurement position to determine the strearawise 
integral length scale of the local turbulence. 

Crosscorrelations generally contained a large positive peak followed by 
a large negative peak. This portion of the signal was approximately anti- 
symmetric about the delay time at zero amplitude. Maximum negative slope, 
corresponding (reference 13) to maximum surface -radiated noise source - 
strength, occurred approximately at this delay time. Normalized cross - 
correlation was taken as the average maximum amplitude of the two peaks. 
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divided by the square root of the product of the maximum values of the auto- 
correlation. These latter maxima occurred at zero delay time. 


Presentation and Interpretation of Data 
Velocity and Turbulence Profiles 

Mean axial velocity profiles half a diameter upstream and downstream of 
the airfoil trailing edges are plotted in figure 46. Flow around the slightly 
curved airfoil surface caused the largest local velocity at these positions 
to exceed the nozzle exhaust velocity. The upper surface shear layer became 
thicker with increased downstream distance. Axial turbulence profiles are 
plotted in figure 4?. Minimum turbulence within the exhaust jet, and maximum 
mean velocity, occurred approximately half a nozzle diameter above the air- 
foil surface. Peak axial turbulence level in the airfoil boundary layer was 
less than 11 percent and occurred where the mean velocity was about 85 per- 
cent of exhaust velocity. Peak axial turbulence level in the lower shear 
layer was about 15 percent. It occurred at the vertical position of the 
trailing edge, where mean velocity was about 70 percent of exhaust velocity. 
Maximum turbulence levels in the upper shear layer were about l4 and 15 per- 
cent for positions upstream and downstream of the trailing edge. These maxi- 
mum levels occurred at a height roughly one nozzle diameter above the trail- 
ing edge, at local mean velocities 65 to 70 percent of nozzle exhaust veloc- 
ity. Peak turbulence levels of 15 percent in the upper and lower shear 
layers at these positions are slightly less than the 17 percent levels shown 
in figure 7 of reference 39 for a position just downstream of the trailing 
edge, for a model with circular exhaust nozzle. 

Crosscorrelation Traces 


Crosscorrelations have been measured between static pressures on the 
upper surface of this USB model with far -field acoustic pressures below the 
model. Those crosscorrelation curves, shown in figure 7b of reference l8, 
had negative peaks followed by positive peaks . From the analysis given in 
reference 13, surface pressures on one side of an airfoil crosscorrelated 
with dipole acoustic pressures in the far field on the same side of that air- 
foil should cause a positive peak followed by a negative peak. Maximum nega- 
tive slope and zero amplitude occur at a delay time which, for a noise source 
location, is equal to the acoustic travel time r/a. The change to a USB con- 
figuration, with acoustic pressures beneath the airfoil occurring l 80 ° out of 
phase with those above the airfoil, produces one reversal of sign. However, 
in isentropic flow an increase of static pressure corresponds to a decrease 
of flow velocity. Use of a hot wire rather than a pressure transducer 
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therefore reverses the signs. Thus a crosscorrelation between a hot wire 
above an airfoil and within a flow that causes dipole noise, and a far field 
microphone below the airfoil, should produce' a positive peak followed by a 
negative peak. In contrast, directly radiated or refracted quadrupole noise 
should produce a single negative peak centered at delay time r/a (references 
4o and 4l) . 

Crosscorrelation signals are plotted in figure 48 for hot-wire locations 
in the airfoil boundary layer and lower shear layer. Amplitudes of the orig- 
inal data traces were replotted such that all curves have the same scale for 
normalized crosscorrelation coefficient. All four curves have the same basic 
shape of a positive peak followed by a negative peak, with maximum negative 
slope at zero amplitude, corresponding to dipole noise. For the two posi- 
tions in the airfoil boundary layer, the zero crossing with maximum negative 
slope occurred at times which exceeded the time required for an acoustic wave 
to travel from the model to the far -field microphone. The zero-crossing time 
approached this acoustic travel time as the trailing edge was approached. 
Similarly, the zero-crossing time measured half a diameter downstream of the 
trailing edge was smaller than the acoustic travel time. This behavior would 
occur if the hot wires were detecting turbulence which maintains its identity 
as it is convected along the airfoil and into the near wake, but causes sur- 
face-radiated noise only as it passes close to the trailing edge. A convec- 
tion velocity can bo inferred from the change of zero-crossing times with 
streamwise distance. This convection velocity is about 60 percent of nozzle 
exhaust velocity, and is approximately equal to the mean velocity at these 
locations of maximum turbulence (figure 46). Approximately this same con- 
vection velocity ratio was reported in reference 40 for tests of a USB con- 
figuration having a 10:1 slot nozzle rather than the circular nozzle used 
with the tests reported herein. 

The crosscorrelation trace which was measured half a diameter downstream 
of the trailing edge has a negative peak at approximately the far-field 
acoustic travel time, preceded by a positive peak. The negative peak corre- 
sponds to quadrupole noise radiated directly to the far field by the local 
high-intensity turbulent flow. As previously mentioned, the combination of a 
positive and negative peak at this position corresponds to surface -radiated 
noise generated by a portion of this turbulence at an earlier time. The 
measurement at one diameter downstream has a positive peak followed by a 
negative peak, followed in turn by a positive peak at the acoustic delay time. 
Here the first pair of peak levels occurs at a time consistent with genera- 
tion of surface -radiated noise near the trailing edge, followed by convection 
of turbulence within the lower shear region to the hot wire. Interpretation of 
the positive peak at the acoustic travel time is not obvious. The analysis 
given in reference 4o showed that a positive second derivative of the 
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crosscorrelation corresponded to direct acoustic radiation from the turbu- 
lence-measuring hot wire to the far -field microphone. Maximum positive sec- 
ond derivative corresponds to zero slope and minimum amplitude. However, 
positive rather than negative peaks of the crosscorrelation were reported in 
reference 4l as representing directly radiated quadrupole noise. Perhaps 
the positive peak at the far -field acoustic travel time then represents 
quadrupole noise from the lower shear layer, and the preceding positive and 
negative peaks represent surface -radiated noise generated when the measured 
turbulence had convected past the trailing edge. 

The largest peak values for normalized crosscorrelation in the lower 
shear layer were those for the measurement location half a diameter upstream 
of the leading edge. However, a comparison of these peak amplitudes alone 
does not give a direct measure of the relative strength of surface -radiated 
and volume -radiated noise sources (references 13, 4o, and 4l) . 

Crosscorrelations between the far-field acoustic pressure and 
velocities in the upper shear region are given in figure 49. These traces 
contain a weak positive peak followed by a strong negative and then a strong 
positive peak. If the two largest peaks are considered to be the dominant 
feature, they represent the downstream convection of a flow disturbance of 
opposite sign to that which produced the far-field noise. The delay times 
at the zero crossing with positive slope are about a millisecond larger than 
those for zero crossings shown in figure 48. If these flow disturbances lie 
in planes skewed about 45° to the jet axis, with their slower-moving portion 
in the upper shear layer, the observed differences in delay times would 
occur . 


Another viewpoint might be that the first positive peak followed by the 
large negative peak is the important feature. These portions of the cross- 
correlations have zero crossing times approximately equal to those in figure 
48 for the airfoil boundary layer and the lower shear region. The large- 
scale structure would then be approximately perpendicular to the exhaust jet 
as with an isolated jet. Interpretation of the final positive peak, and in 
particular its decrease of occurrence time with increasing downstream dis- 
tance, then becomes difficult. The most likely interpretation of these data 
for the upper shear layer is that the upper shear layer of an axisymmetric- 
nozz] e USB is related to noise radiation in a complicated manner that cannot 
be readily approximated as a sum of distinct volume -radiated and surface- 
radiated components . 
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Results of Crosscorrelation 


Spatial distributions of axial integral length scale v/ithin the surveyed 
region are plotted in figure 50a. This length scale generally was about l/3 
the nozzle diameter. It increased to about 45 percent of this diameter in 
the upper shear layer downstream of the trailing edge. Thus the positions of 
maximum turbulence intensity are associated with large-scale structure of the 
Jet turbulence similar to that of an isolated jet. 

Maximum normalized crosscorrelation coefficients between the jet veloc- 
ity and far-field acoustic pressure are shown in figure 50b. They were evalu- 
ated for the largest adjacent positive and negative crosscorrelation peaks in 
figures 48 and 49. Maximum normalized crosscorrelation coefficient ranged from 
0.08 to 0.10 for most locations. The only significantly larger value, 0.13, 
was measured in the airfoil boundary layer half a diameter upstream of the 
trailing edge. Normalized crosscorrelations of upper surface pressures and 
far field acoustic pressures for this USB configuration were shown in figure 
lOd of reference 18 for a somewhat higher exhaust velocity (l25 m/sec). The 
measured value of nearly 0.4 at a position nearly half a diameter upstream of 
the trailing edge is much larger than the largest value measured with the hot 
wire. Evidently, the hot wire measured far more details of the local 
unsteady flow than the area average provided by a surface pressure trans- 
ducer . 

From the crosscorrelation amplitudes, and the delay times for zero 
crossings with large -amplitude slopes, it is concluded that USB configura- 
tions with axisymmetric nozzles and short flap lengths have two important 
noise processes. One is convection of airfoil upper-surface boundary layer 
turbulence past the neighborhood of the trailing edge, generating a surface - 
radiated edge noise. The other is direct radiation of quadrupole noise from 
the high-turbulence lower shear layer at moderate distances downstream of the 
trailing edge. 

It should be noted that two concurrent studies (references 4l and 42) 
have investigated noise radiation from USB slot nozzle configurations. Both 
experimental programs used an aspect ratio 10 slot nozzle and a flat -plate 
wing having zero flap deflection. Both wing models were longer than the jet 
potential core, unlike the situation for the configuration described herein. 
Both programs used crosscorrelations between flow-field hot wires and far- 
field microphones. In reference 4l it was concluded that the dominant process 
was quadrupole noise radiated directly from a small region of high turbulence, 
several slot heights downstream of the trailing edge. However, for the 
nearly identical configuration and test program described in reference 42, 
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trailing-edge noise caused by turbulence in the upper shear layer was 
reported to dominate. If noise radiation from USB configurations can be 
approximated by a sum of several simplified noise components, associated 
with discrete source locations, the location of the dominant source cannot 
be conclusively identified for all cases. 
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UNDER-THE- WING NOISE MECHANISM 

(UTW) CONFIGURATION 


UPPER- SURFACE -BLOWING 
lUSB) CONFIGURATION 


VANE 



FLUCTUATING 
LIFT NOISE, 
INFLOW NOISE 





FIGURE 1. - SKETCH OF DIRECTIVITIES FOR DIFFERENT MECHANISMS 
OF EXTERNALLY BLOWN FLAP NOISE 



ONE-THIRD-OCTAVE SOUND PRESSURE LEVEL, dB 




ONE-THIRD- OCTAVE BAND CENTER FREQUENCY, Hz 


CIRCLE SYMBOLS ARE SUM OF MEASURED SOUND ON SIDE SHEILED FROM JET, PLUS 
MEASURED NOZZLE-ALONE SOUND INCREASED 3dB BY REFLECTION FROM WING 


FIGURE 2. - COMPARISON OF FAR-FIELD SPECTRA IN THE SCRUBBED AND UNSCRUBBED 
DIRECTIONS FROM AN UNDEFLECTED EXTERNALLY BLOWN WING AT 250 M/SEC VELOCITY 
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(b.) TRIPLE SLOTTED FLAP MODEL APPROACH CONFIGURATION 


(c.) 




SLOTLESS WING WITH TRIPLE SLOTTED FLAP APPROACH CONTOUR 


FIGURE 3. -UNDER-THE-WING CONFIGURATIONS FOR 
EVALUATION OF EBF NOISE PREDICTION METHODS. 
MODEL CONFIGURATIONS ARE DRAWN TO SCALE. 
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FIGURE 4. -UPPER SURFACE BLOWING CONFIGURATIONS FOR 
EVALUATION OF EBF NOISE PREDICTION METHODS. 
MODEL CONFIGURATIONS DRAWN TO SCALE. 
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OVERALL SOUND PRESSURE LEVEL, dB OVERALL SOUND PRESSURE LEVEL, dB 




POLAR ANGLE, e DEG 


FIGURE 5. - CALCULATED AND MEASURED OASPL DIRECTIVITY IN THE FLYOVER PLANE 

FOR UTW MODEL APPROACH CONFIGURATION 
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OVERALL SOUND PRESSURE LEVEL AT SIDELINE ANGLE MINUS OVERALL SOUND 
PRESSURE LEVEL IN F LYOVER PLANE, (OASPL}^ _ (OASPL)^ = q, dB 
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(a) 60° sideline ANGLE 
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POLAR ANGLE, e, DEG 
(b) 85° sideline angle 

FIGURE 6. - EFFECT OF SIDELINE ANGLE ON CALCULATED 
AND MEASURED OASPL DIRECTIVITY FOR UTW MODEL 
APPROACH CONFIGURATION 

8i 


1/3 OCTAVE SPL-OASPL clB 1/3 OCTAVE SPL-OASPL, rIB 



FLYOVER PLANE 



STROUHAL NUMBER, fD/V 
1b) 85° SIDELINE PLANE 


FIGURE 7. - EFFECT OF POLAR ANGLE AND EXHAUST VELOCITY 
ON CALCULATED AND MEASURED NORMALIZED 1/3 OCTAVE 
SPECTRA FOR UTW MODE L APPROACH CONFIGURATION 
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OVERALL SOUND PRESSURE LEVEL, dB OVERALL SOUND PRESSURE LEVEL, dB 




FIGURE 8. - CALCULATED AND MEASURED OASPL DIRECTIVITY IN THE 
FLYOVER PLANE FOR UTW MODEL TAKEOFF CONFIGURATION 





OVERALL SOUND PRESSURE LEVEL AT SIDELINE ANGLE MINUS OVERALL SOUND PRESSURE 
LEVEL IN FLYOVER PLANE, (OASPL)0-(OASPL)0 = q, dB 




FIGURE 9. - EFFECT OF SIDELINE ANGLE ON CALCULATED AND MEASURED OASP 
DIRECTIVITY FOR UTW MODEL TAKEOFF CONFIGURATION 
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IGURE 10. - EFFECT OF POLAR ANGLE AND EXHAUST VELOCITY ON CALCULATED AND 
MEASURED NORMALIZED 1/3 OCTAVE SPECTRA FOR UTW MODEL 
TAKEOFF CONFIGURATION 
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1/3 OCTAVE SOUND PRESSURE LEVEL, dB 



FIGURE 12 COMPARISON OF MEASURED SPECTRA WITH THOSE CALCULATED 
BY UTRC METHOD FOR POSITION DIRECTLY BENEATH 
SLOTLESS WING 
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FIGURE 13. - EFFECT OF POLAR ANGLE AND EXHAUST 

VELOCITY ON CALCULATED AND MEASURED NORMALIZED 
1/3 OCTAVE SPECTRA FOR UTW MODEL SLOTLESS WING 
IN FLYOVER PLANE 
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iURE 14.- CALCULATED AND MEASURED OASPL DIRECTIVITY IN THE FLYOVER PLANE 
FOR QCSEE USB 1/11.5 SCALE MODEL, TAKEOFF CONFIGURATION 
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FIGURE 15. - EFFECT OF POLAR ANGLE ON CALCULATED AND MEASURED SIDELINE 
DIRECTIVITY FOR QCSEE USB 1/11.5 SCALE MODEL, TAKEOFF CONFIGURATION, 218 M/SE 
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(a) FLYOVER PLANE 



STROUHAL NUMBER ,fD/V 
(b) 63° SIDELINE PLANE 


URE 16. - EFFECT OF POLAR ANGLE AND EXHAUST VELOCITY ON CALCULATED AND 
1EASURED NORMALIZED 1/3 OCTAVE SPECTRA FOR QCSEE 1/11.5 SCALE USB MODEL 

TAKEOFF CONFIGURATION 
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FIGURE 17. - EFFECT OF FLAP LENGTH ON DIRECTIVITY AND NORMALIZED SPECTRJ 
OF NOISE FROM QCSEE USB 1/11.5 SCALE MODEL, TAKEOFF CONFIGURATION, 218 M/sJ 

EXHAUST VELOCITY 1 
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POLAR ANGLE, d, DEG 

URE 18. - CALCULATED AND MEASURED OASPL DIRECTIVITY IN THE FLYOVER PLANE FOR 
QCSEE USB 1/11.5 SCALE MODEL APPROACH CONFIGURATION 
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FIGURE 19. - EFFECT OF POLAR ANGLE ON CALCULATED AND MEASURED SIDELINE 
DIRECTIVITY FOR QCSEE USB 1/11.5 SCALE MODEL, APPROACH CONFIGURATION, 190 M/SI 

EXHAUST VELOCITY 
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(b) 630 SIDELINE PLANE 


sURE 20.- EFFECT OF POLAR ANGLE AND EXHAUST VELOCITY ON CALCULATED AND 
ASURED NORMALIZED 1/3 OCTAVE SPECTRA FOR QCSEE 1/11.5 SCALE USB MODEL 

APPROACH CONFIGURATION 
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FIGURE 21. - CALCULATED AND MEASURED OASPL DIRECTIVITY 
IN THE FLYOVER PLANE FOR TF-34 USB 1/18.5 SCALE MODEL 
WITH SLOT NOZZLE AND 40° DEFLECTION SHORT FLAP 
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JRE 22. - EFFECT OF POLAR ANGLE AT TWO EXHAUST VELOCITIES ON CALCULATED 
AND MEASURED NORMALIZED 1/3 OCTAVE SPECTRA FOR TF-34 USB 1/18.5 
SCALE MODEL WITH SLOT NOZZLE AND 40° DEFLECTION SHORT FLAP 
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OVERALL SOUND PRESSURE LEVEL, dB 




FIGURE 23.- CALCULATED AND MEASURED OASPL DIRECTIVITY IN THE 
FLYOVER PLANE FOR USB WITH CIRCULAR NOZZLE AND VANE 
DEFLECTOR IN TAKEOFF CONFIGURATION 
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FIGURE 24.- EFFECT OF POLAR AMGLE ON CALCULATED AND MEASURED 
SIDELINE DIRECTIVITY FOR USB WITH CIRCULAR NOZZLE AND VANE 
DEFLECTOR IN TAKEOFF CONFIGURATION 
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FIGURE 25._caLCULATED AND MEASURED OASPL DIRECTIVITY 
IN THE FLYOVER PLANE FOR USB WITH CIRCULAR NOZZLE 
AND VANE DEFLECTOR IN APPROACH CONFIGURATION 
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FIGURE 26. -EFFECT OF POLAR ANGLE AND EXHAUST VELOCITY ON CALCULATED AND 
MEASURED NORMALIZED 1/3 OCTAVE SPECTRA IN FLYOVER PLANE FOR USB WITH 
RCULAR NOZZLE AND VANE DEFLECTOR IN TAKEOFF AND APPROACH CONFIGURATIONS 
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(a) 60° SIDELINE ANGLE 



STROUHAL NUMBER, fD/V 


(b) 85° SIDELINE ANGLE 


FIGURE 27, -EFFECT OF EXHAUST VELOCITY AND SIDELINE ANGLE ON CALCULATED 
AND MEASURED NORMALIZED 1/3 OCTAVE SPECTRA AT 90° POLAR ANGLE FOR USB 
WITH CIRCULAR NOZZLE AND VANE DEFLECTOR IN TAKEOFF CONFIGURATION 


102 





OVERALL SOUND PRESSURE LEVEL, dB OVERALL SOUND PRESSURE LEVEL, dB 


130 


120h 


iioh 


WING WITH 
PLUG FAIRINGS 



UTRC 




NASA LEWIS 
CALCULATED 


EXHAUST VELOCITY 
255 M/ SEC 


100 1 



40 80 

POLAR ANGLE, 6 , dB 


120 


160 



FIGURE 28. - CALCULATED AND MEASURED OASPL DIRECTIVITIES IN 

THE FLYOVER PLANE FOR LARGE DOUBLE SLOTTED WING 
WITH PLUG FAIRINGS IN SLOTS 
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FIGURE 29.- CALCULATED AND MEASURED OASPL DIRECTIVITY IN 
FLYOVER PLANE FOR UTW MODEL SLOTLESS WING 
HAVING CONTOUR OF DOUBLE SLOTTED FLAP 
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FIGURE 30. -CALCULATED AND MEASURED OASPL DIRECTIVITY IN 
FLYOVER PLANE FOR 45° DEFLECTION UTW SLOTLESS 
WING AT THREE FLAP LENGTHS AND 270 M/SEC EXHAUST 
VELOCITY 
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(b) APPROACH FLAP CONFIGURATION 

FIGURE 31, - COMPARISON OF CALCULATED AND MEASURED OASPL 
DIRECTIVITY FOR LARGE-SCALE UTW MIXER NOZZLE 
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(b) APPROACH FLAP DEFLECTION 


SURE 32. - COMPARISON OF CALCULATED AND MEASURED OASPL DIRECTIVITY FOR 
ENGINE-IN-FRONT-OF-WING CONFIGURATION DATA DECREASED 3dB 
FOR COMPARISON WITH FREE-FIELD PREDICTIONS 
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(a) TAKEOFF FLAP DEFLECTION, (b) APPROACH FLAP DEFLECTION, 

100° MICROPHONE ANGLE 80° MICROPHONE ANGLE 


FIGURE 33. - COMPARISON OF CALCULATED AND MEASURED SPECTRA FOR 
ENGINE-IN-FRONT-OF-WING CONFIGURATION. DATA DECREASED 3 dB 
FOR COMPARISON WITH FREE-FIELD PREDICTIONS 
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CURE 34. - SKETCH OF FULL SCALE QCSEE EIMGIIME EXHAUST NOZZLE AND EXTERNALLY 
BLOWN FLAP CONFIGURATIONS 
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URE 36. - CALCULATED FREE FIELD DIRECTIVITY OF FULL SCALE QCSEE UTW APPROACH 
CONFIGURATION AT 100M RADIUS AND 161 M/SEC EXHAUST VELOCITY BY 
UTRC PREDICTION METHOD 


111 




1/3 OCTAVE SOUND PRESSURE LEVEL,dB 1/3 OCTAVE SOUND PRESSURE LEVEL, dB 



50 100 200 500 1000 2000 5000 10000 

1/3 OCTAVE BAND CENTER FREQUENCY.Hz 



50 100 200 500 1000 2000 5000 10000 


1/3 OCTAVE BAND CENTER FREQUENCY.Hz 


FIGURE 37. -CALCULATED 1/3 OCTAVE SPECTRA AT 90° POLAR ANGLE FOR FULL SCALE 
QCSEE UTW IN TAKEOFF AND APPROACH CONFIGURATIONS AT100M RADIU 
BY UTRC PREDICTION METHOD 
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POLAR ANGLE, 0, DEG 

CALCULATED FREE FIELD DIRECTIVITY OF FULL SCALE QCSEE USB APPROACH 

CONFIGURATION AT 100 M RADIUS AND 190 M/SEC EXHAUST VELOCITY 
BY UTRC PREDICTION METHOD 



EXHAUST VELOCITY EXPONENT, n, 

WHERE OASPL= 10 n LOG Vg + K OVERALL SOUND PRESSURE LEVEL, OASPL, dB 




FIGURE 39. - CALCULATED FREE FIELD DIRECTIVITY OF FULL SCALE QCSEE USB TAKE 
CONFIGURATION AT 100 M RADIUS AND 220 M/SEC EXHAUST VELOCITY 
BY UTRC PREDICTION METHOD 
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FIGURE 40.- CALCULATED 1/3 OCTAVE SPECTRA AT 90° POLAR ANGLE FOR 
FULL SCALE QCSEE USB IN TAKEOFF AND APPROACH 
CONFIGURATIONS AT 100 M RADIUS BY UTRC PREDICTION 
METHOD 
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FIGURE 41. -CALCULATED EFFECT OF 41 M/SEC (80 KNOT) FLIGHT VELOCITY Ol| 
FLYOVER SPECTRUM AT OVERHEAD POSITION AND 100 M ALTITUDE FOR QCSEeI 
USB CONFIGURATIONS BY UTRC PREDICTION METHOD I 
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CURE 42. -CALCULATED EFFECT OF 41 M/SEC (80 KNOT) FLIGHT VELOCITY ON 
FLYOVER SPECTRUM AT OVERHEAD POSITION AND 100 M ALTITUDE FOR 
QCSEE UTW CONFIGURATIONS BY UTRC PREDICTION METHOD 
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ENGINE IN FRONT OF THE WING, CONFIG = 3 


FIGURE 43. - GEOMETRIC INPUT FOR EXTERNALLY BLOWN FLAP CONFIGURATIONS 





FIGURE 44. -FLOW CHART FOR EBF NOISE PREDICTION COMPUTER PROGRAM 
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FIGURE 45. 


SKETCH OF EXHAUST JET FLOW REGIONS FOR UPPER 
SURFACE BLOWING EXTERNALLY BLOWN FLAP 


RATIO OF MEAN AXIAL VELOCITY TO NOZZLE VELOCITY, U/Uj 



DISTANCE NORMAL TO TRAILING EDGE POSITION, Y/D 


IGURE 46. - MEAN VELOCITY PROFILES UPSTREAM AND DOWNSTREAM OF UPPER 
SURFACE BLOWING MODEL TRAILING EDGE. SPANWISE CENTERLINE, 
100 M/SEC NOZZLE VELOCITY 
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RMS AXIAL TURBULENCE LEVEL, U/U 



DISTANCE NORMAL TO TRAILING EDGE POSITION, Y/D 


FIGURE 47- AXIAL TURBULENCE VELOCITY PROFILES UPSTREAM AND DOWNSTREAM 
OF UPPER SURFACE BLOWING MODEL TRAILING EDGE, SPANWISE 
CENTERLINE, 100 M/SEC NOZZLE VELOCITY 
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FIGURE 49. -CROSSCORRELATIOIMS BETl 
AND ACOUSTIC FAR FIELD 
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(b) MAXIMUM NORMALIZED CROSSCORRELATION COEFFICIENTS BETWEEN 
EXHAUST JET VELOCITY FLUCTUATIONS AND FAR FIELD NOISE 


SURE 50. -SPATIAL DISTRIBUTION OF INTEGRAL SCALE LENGTH. AND CORRELATION 
BETWEEN LOCAL VELOCITY AND FAR FIELD NOISE, NEAR UPPER SURFACE 
BLOWING MODEL TRAILING EDGE. SPANWISE CENTERLINE, 100 M/SEC 
NOZZLE VELOCITY 
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